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THE FIVE-DIMENSIONAL GEOMETRY OF THE 
CURVATURE TENSOR IN A RIEMANNIAN J, 


By H. 8. RUSE (Southampton) 
[Received 27 February 1944] 


1. Introduction 
Ir X' is a contravariant vector at a given non-singular point P, (2*), 
of a Riemannian V, of fundamental tensor g;;, then the X* may be 
interpreted as the homogeneous coordinates of a point in the pro- 
jective (n—1)-space S,_, at infinity in the centred affine tangent- 
space 7), at P. If 
pi = Xtyi—xiyé 

are Pliicker coordinates of the line joining points X‘, Y‘ of S,_, (the 
p" therefore being the components of a simple bivector of V, at P), 
then the lines for which 

= 0, (1.1) 
where R;;,, is the covariant Riemann tensor (skew in ¢, j and in k, 1), 
form a quadratic complex in S,_,, the Riemann complex. 

In a previous paper (Ruse, 6) it was shown that, for n = 3, the 
theory of the Riemann complex reduces to a simple one of conics 
in a projective S,, and an introductory discussion was given of the 
case n = 4. It is the purpose of this paper to continue that discus- 
sion. Hereafter all Latin indices will run from 1 to 4. 

For n = 4 the theory is one of a quadratic complex in a projective 
S;, together with a non-degenerate quadric, the fundamental quadric, 
of point-equation g,; X‘X’ = 0, of tangential equation g“u;u; = 0, 
and of line-equation 

Jinn = 0, or = 0, (1.2) 

and where °p;; are the dual coordinates of the line p”, defined by 

= P™. (1.4) 

Here the } is a factor of proportionality chosen for convenience, and 

€:j« 18 the dualizing tensor, or dualizor, of components +g, 0. The 

dualizors €;;,,, €™ (= +1/vg, 0) are imaginary if P is a point of V, 

at which g < 0. When the fundamental quadric is taken in either 

of the forms (1.2), it is to be regarded as defined by the quadratic 
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complex of lines touching it, (1.2) being dual equations of this special 


complex because ogitl — Jetimnglipag 


(1.5) 
by Jacobi’s theorem applied to the second-order minors of the adju- 


gate of the determinant |g;;|. 
The dual equation of the Riemann complex is 


= 0, (1.6) 


RPijkl — 1 ijmn kl 
where = (1.7) 


Raising and lowering suffixes by means of the fundamental tensor 
corresponds in S, to taking the polar with respect to the fundamental 
quadric. Thus, if € is a point of S,, then 

= 9ij gi 
is its polar plane, of equation €; X* = 0 in current point-coordinates 
X‘. For this reason the operation of raising and lowering suffixes 
by means of the fundamental tensor will sometimes be referred to 
briefly as polarizing; and the operation of raising and lowering pairs 
of skew suffixes by means of the dualizors, as in (1.4), (1.5), (1.7), 
will be called dualizing. The two operations are geometrically equi- 
valent, and also algebraically equivalent except perhaps for a multi- 
plicative scalar, when the geometrical configuration represented by 
the operand is self-polar with respect to the fundamental quadric. 
Thus, for example, the line p“” is self-polar when, and only when, 
there exists a scalar p such that 

= PPij- 
When p = 1, the operations are algebraically equivalent. It is to be 
observed that, in all cases, the two operations are commutative 
because—cf. (1.5)— 

The word polar by itself will always refer to the polar with respect 
to the fundamental quadric. 

The contravariant Riemann tensor R“* and its dual °R;,;,,; define 
the quadratic complex, of dual equations 

= 9, = 0, (1.9) 


polar to the Riemann complex. °R;,,, and R® are obtainable from 
one another by dualizing, and are respectively obtained from °R* 
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and R;,,, by polarizing. The unsatisfactoriness of this notation from 
the geometrical point of view was commented upon in §5 of the 
previous paper quoted above. 

The Riemann complex is self-polar when 


— Risk, (1.10) 


It was shown in the previous paper that the scalar « must have one 
of the values +1. When « = +1, (1.10) is equivalent to the Einstein 


condition 


2. Ennuplet coordinates in S, 

Let hi, = (hi,...,h4) be the components of an orthogonal ennuple 
at P in V,. The letters a, 6,..., g will always be used for ennuplet 
suffixes and #, j, k,... for proper tensor suffixes. We suppose that the 
h’s satisfy the usual relations, and in particular 


= gi. (2.1) 


By taking a simple summation on the right-hand side, omitting the 
‘indicators’ e, (= +1) which correspond to the signature of the 
fundamental quadratic form (Eisenhart, 2, ch. IIT), we are treating 
the hi, as though the fundamental form were positive-definite at (2*). 
But the possibility of an indefinite metric need not be excluded, 
since it is of no consequence in the present theory if some of the hi 
are imaginary. 

Equation (2.1)-may be written 

hi, = g", (2.2) 

where g® is equal to the Kronecker symbol 5”. [This notation has 
a number of advantages: for example, the possibility of an indefinite 
metric may be provided for by taking the h’s real and g® = e,8% 
(not summed) instead of assuming, as we are doing, that some of 
the h’s may be imaginary. | 

If h = det|hi|, it follows from (2.2) that h? = 1/g, and we suppose 
that the signs of the h’s are so chosen that 


h>0. (2.3) 
If g,» is equal to the Kronecker symbol 6,,, then also 
Gar hi = 
where [h?] is the matrix reciprocal to [hi]. Also 
he = 9% gi; hj. 


| 
Ri; = Bi; 
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Thus g” and g,,, may be used to raise and lower ennuplet suffixes. 
The transformation = 
from tensor to ennuplet components of a contravariant vector corre- 
sponds in the projective space S, to a change of homogeneous 
coordinates, the vertices of the new tetrahedron of reference being 
the points whose coordinates in the X‘-system are hj,..., hj. The X¢ 
are scalars for transformations of the coordinates x of V,. Ennuplet 
components of tensors are distinguished from coordinate-components 
merely by having suffixes a, b,... instead of ¢, j,..., the same central 
letter being used in each case. The ennuplet components of the 
fundamental tensor itself are g;;hi,hj, and this is equal to g,,. The 
equation of the fundamental quadric in S, in ennuple. coordinates 


is thus P 
Ja A°X® = 0, ie. > (X* = 0. 
a=1 
The new tetrahedron of reference is therefore self-polar with respect 
to the fundamental quadric. 
The ennuplet components of the Riemann tensor are 
Ravea = hi, hj hk hi, Raved = Riikihe hé 


The ennuplet dualizors ¢€,,,q and <«“’°? both have components +1, 0, 


abed 


as follows from (2.3). 


3. Five-dimensional representation 

The representation of the lines of a projective 3-space by the 
points of a 4-quadric in a projective 5-space is too well known to 
need exposition (see, e.g., Sommerville, 8, 343), but it is necessary 
to describe it briefly in order to explain the notation used in the 
present paper. 

Let g@ (= —q*) be the ennuplet components of a linear complex 
in S,, its equation in current line-coordinates °p,,, being g@°p,, = 0, 
or dually, °¢,,p* = 0. Write 

(3.1) 
or, in one equation, 
(a = I, 2...., 6). (3.2) 
The coefficients y*, are those of the linear transformation (3.1). 
Thus, for example, 
Ys = —Ys2 = (1,0, 0, 0, 0, 0). 
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Greek suffixes will always run from 1 to 6. The g* may be inter- 
preted as the homogeneous coordinates of a point in a projective 
space S;. The transformation (3.2) may also be written 


q” (3.3) 
where = = (1,0, 0,0, 0,0), etc. 
As it stands at present the transformation amounts merely to the 


replacement of pairs of skew suffixes by single suffixes according to the 
scheme 


4 
6 


By (3.2) and (3.3) it quickly follows that 
bye = 88, (3.5) 
and that fr yz, = 52, (3.6) 
where 5% is the generalized Kronecker symbol 8233—898°. Write 
€ap = Ver (3.7) 
Then €,g, e** are symmetric in a, B, and 
= 53. (3.8) 
Also Eabed = Yar Yea (3.9) 
Equations (3.7) are a disguised form of the statement—cf. (3.4)—that 


ab = —ba= 23 31 12 14 24 7 


€14 = €41 = 9314 = 1, 
= €42 = €3124 = 1, 
€34 = €43 = €1934 = 1, 


with €4g = 0 otherwise, and similarly for «*8. Thus 


= (a row, column), (3.10) 


I and O being the unit and null 3-by-3 matrices respectively. If 
Gap = 9° = (3.11) 

where, as usual, Javea = Jac Ioa—Gaa 

then 9°"9yg = (3.12) 

Because g,», g@ are both equal to the Kronecker delta, it follows that 


10 
Gap = Fup, gr = ive. 


(3.13) 
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The formalism may seem a trifle elaborate for so simple a trans- 
formation as that implied by (3.4). The reason for its introduction 
is twofold: it makes possible the use of Greek suffixes as applying 
to S;, and also frees the theory from the restriction to the particular 
coordinate-systems implied by the choice of the ennuple hi, in S, and 
by the transformation (3.1). Thus, if we write 

Vii = hj, 
then, for fixed «, y#; is a covariant six-vector of V, which defines 
a linear complex in S3, and (3.2) may be written 

= 3759". (3.14) 
Moreover, in S; the coordinates g* may be subjected to linear trans- 


formations, say 


and (3.14) then becomes 

q* = 
where vi = 
In other words, the coordinates of S, and 8S; may be transformed 
independently of one another if y%; is treated as a tensor of the type 
indicated by its suffixes. In all proper tensor formulae, such as (3.5), 
(3.6), (3.7), and (3.8), ennuplet suffixes may be replaced by proper 
tensor suffixes, or may be taken to refer to an ennuple other than 
hi, and the Greek suffixes may be taken to refer to an arbitrary 
coordinate-system in S;. For general coordinate-systems the numeri- 
cal values of €,g, eB, Jap» g*? given by (3.10) and (3.13) do not hold. 
When they do hold, that is, when the transformation from S, to S; 
is given by (3.1), the coordinate-systems q* and X®@ in S, and S; 
respectively will be called basic, and the orthogonal ennuple hi, will 
be called the basis of the transformation. 

To pass from a tensor equation in V, or 8, to the corresponding 
equation for S,, it is therefore necessary only to replace pairs of 
skew suffixes (tensor or ennuplet) by single Greek suffixes, and, when 
necessary, to calculate factors of proportionality by reference to the 
basic coordinate-system. 

Now the coordinates p” of a line in S, satisfy the identity 


which, multiplied by 2vg, may be written 
p™ = 0, or, dually, = 0. (3.15) 
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The corresponding equations for S; are 


= 0 (3.16) 
and = 0, (3.17) 
where Pa = PP. (3.18) 


Of these (3.16) is the point-equation in current coordinates p* of the 
non-degenerate 4-quadric in S; the points of which correspond to 
the lines of S8;. Because of (3.8) it is evident that (3.17) is the 
tangential equation of the same 4-quadric in current hyperplane- 
coordinates °p,. For any point p*%, (3.18) defines its polar 4-plane 
with respect to the 4-quadric €,g, which will be called the e-quadric. 
The lines of the special quadratic complex 

= 0 (3.19) 
in S, correspond in S; to the points of the 3-space of intersection of 
the non-degenerate 4-quadric 

Japp"p? = 0 (3.20) 

with the e-quadric. We call this the g-quadric, the term ‘funda- 
mental’ being inappropriate, even though (3.20) is obtained from 
the line-equation (3.19) of the fundamental quadric in S;, because 
of the even more fundamental nature of the e-quadric in S;. By 
(3.12) the tangential equation of the g-quadric is 

= 0. 
By (1.5), 5 (3.21) 
The left-hand side represents the polar of the g-quadric with respect 
to the e-quadric, and the right-hand side is the tangential form of 
the g-quadric. Hence the g-quadric is self-polar with respect to the 
e-quadric. Similarly, since 

(this being merely the statement that the e-tensors are obtainable 
from one another by polarizing), we have, in S,, 

= (3.22) 
Therefore the e-quadric is self-polar with respect to the g-quadric. 
Equations (3.21) and (3.22) are both deducible from 


yp = yp 
—cf. (1.8)—which states that the operations of taking the polars 
with respect to the two quadrics are commutative. 
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The polars of any geometrical configuration in S; with respect to 
the quadrics will be called its e-polar and g-polar respectively. 
The lines of the Riemann complex 


= 0 
in S, correspond in S; to the points common to the 4-quadric 
R, == @ 
and the e-quadric. Here 
The dual equation of the complex in S; is 
and the corresponding equation for S; is 
°RP°p, “pg = 0, 
where — Rs = Rise, 
In the basic coordinate-system (3.24) is 


the (ennuplet) suffix-pairs on the right running according to the 
scheme (3.4). The symmetric matrix (3.28) may be symbolized by 
Rog=[P 
3.29 
where, of the 3-by-3 matrices P,Q, S,S‘, P and Q are symmetric 
and S‘ is the transpose of S. 
By (3.27) °R°? is the 4-quadric (envelope) in S,; polar to the 
4-quadric R,, with respect to the e-quadric. The point-equation of 


-quadric is 
the 4-quadric is ° == ©, 


where the coefficients °RYg are such that 
« being a scalar (or zero if the 4-quadric is degenerate). In general 
the tensor °RYB is not the same as the original tensor Ry, nor even 
a scalar multiple of it: it is so only when the quadric Ry is self-polar 
with respect to the e-quadric, and in that case, by (3.30), 
°RY = pop (p scalar). (3.31) 


(3.23) 
(3.24) 
(3.25) 
(3.26) 
(3.27) 
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In S, this becomes ° Rime Ras = po, (3.32) 
which, as will be seen in § 5, equation (5.5), is the condition that the 
Riemann complex should be special, that is, consist of the tangents 
to a quadric, or, degenerately, of the tangents to a cone or of the 
lines meeting a proper or degenerate conic. 

Similarly, the tangential equation of the quadric R,¢ is 

R*B°p =0, 
where is such that = 
and R*8 is the same as °R*, or is a scalar multiple of it, only when 
the complex in S; is special. 

As it is desirable as far as possible to avoid a confusion of nota- 
tions, no more use will be made of the tensors ° RX, and R*8, which 
correspond in S, to the quadratic complex conjugate to the Riemann 
complex (Salmon, 7, 51). 

In addition to the quadric Rg and its e-polar °R*’, we have, in 
S,;, the 4-quadric (envelope) defined by 

which corresponds to the polar complex R’* in S,, and the 4-quadric 
(locus) defined by 


which corresponds in S, to the dual °R;;,, of the polar complex R&™, 
The former is the g-polar of the original 4-quadric R,g. The latter 
is the g-polar of the 4-quadric °R°?, which is itself the e-polar of R,g. 
The envelope R*? has a point-equation R*,p*p? = 0, and the locus 
°R,g @ tangential form °R**, but, once again, we omit these from 
consideration. 

To sum up: In S; we have, fundamentally, three 4-quadrics €xp> 
Jag: Rag, the first being the one whose points correspond to the lines 
of S;. The quadrics €,g and g,g are non-degenerate, and each is self- 
polar with respect to the other. R,g may or may not be degenerate 
according to the nature of the quadratic complex in 8,, the lines of 
which are represented in S; by the intersection of ¢,g and R,g. From 
the quadric R,g we get, taking polars, 

(i) its e-polar (envelope) 
(ii) its g-polar (envelope) R*?; 
(iii) the g-polar °R,, of °R*’, which is the same locus as the e-polar 
of R*?, because the operations of taking the e- and g-polars 
are commutative. 


4 
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We shall not need to consider the polars of the e- and g-quadrics 
with respect to the R-quadrics. 

It is perhaps not surprising that such a wealth of geometrical 
material leads to a fairly extensive theory, to which the present 
paper is intended to serve as an introduction. 


4. The Segre characteristics of the Riemann tensor at a 
point of J, 
Let g* be a point of S;. Its polar 4-planes with respect to the 
R- and e-quadrics have coordinates R,,q°, €,gq°, and are the same if 


= 


that is, if (Rag—Aeag)g’ = 0, (4.1) 
A being a (scalar) root of the equation 

det | (4.2) 
of the sixth degree. When the roots are all distinct, there is one 
point g* corresponding to each root, the six points so obtained being 
the vertices of the common self-polar simplex of the two quadrics. 
In S, the equation corresponding to (4.1) is 

= 9, (4.3) 
and q is a principal linear complex of the quadratic complex R;,,,. 
In V, each q” is a contravariant six-vector. 

Whether the roots of (4.1) are distinct or not, the quadratic com- 
plex R,;,; may be classified by means of the Segre characteristic of 
the matrix pencil [R,g—Ae,g| formed from the exponents of the 
elementary divisors. This provides the usual classification of the 
complex as one of 49 different species (Jessop, 3, 230-2; or Zindler, 
10, 1128-31), or as one of 8 degenerate forms (Zindler, ibid., 1133). 
In my previous paper (Ruse, 6, § 4) it was seen that there is no 
corresponding theory for a V3. 

This classification takes no account, however, of the relation of 
the quadratic complex to the fundamental quadric. This relation 
may be similarly symbolized by the Segre characteristic of the pencil 
[2.g—H9 xp] This corresponds in a general sense to the classification 
of a Vz; by means of the relation between the Riemann conic, or 
equivalently of the Ricci conic, and the fundamental conic (Ruse, 6). 
This simultaneous use of the Segre characteristics of | R..g—Ae,g] and 
[Rg—vGop], Which will be called the e- and g-discriminants respec- 


i 
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tively, obviously gives a more detailed classification of the Rieman- 
nian V, than is afforded by the roots of the equation 


appearing in the theory of Ricci principal directions. 

The two characteristics will be called the «- and g-characteristics 
respectively, and the roots of the discriminant equations the «- and 
g-roots. 

Although the characteristics may be used for a fairly detailed 
description of the nature of a V, at a point, it should be noted that 
they do not by themselves give full information about the Riemann 
tensor at that point: they take into account only equalities among 
the e- and g-roots, and the distribution of those roots among the 
minors of the corresponding matrices. Yet the roots may be con- 
nected by relations other than that of simple equality, as, for 
example, in the case of the harmonic complex (see § 5 below), which, 
though a complex of a particular sort, may have any one of a number 
of e- and g-characteristics, including [111111] (Jessop, 3, 359). 

That both characteristics enter naturally into a discussion of the 
properties of the Riemann tensor in a V, has not been fully recog- 
nized, or at any rate not emphasized, by previous writers. Struik (9) 
and Lamson (4) refer explicitly only to the e-characteristic, though 
Lamson, who confines his discussion to the Einstein V, (R,; = kg;;), 
takes full account of the geometrical relationship of the Riemann 
complex to the fundamental quadric. Churchill (1, §11), on the 
other hand, takes explicit notice only of the g-discriminant. 


5. Special cases 

Consider the case when the V, is of constant curvature. Then at 
Rijn = 
Hence—see (1.2)—the Riemann complex coincides with the special 
complex of tangents to the fundamental quadric, and is therefore 
(Jessop, 3, 211) of e-characteristic [(111)(111)]. Now every special 
quadratic complex has the same e-characteristic, so that a space for 


which 
= (5.1) 


where a,; is any symmetric tensor of non-zero determinant, is also 
of e-characteristic [(111)(111)]. But any V, immersible in a flat space 


det | R;;—pg;;| = 0 (4.4) 
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of five dimensions is of this type provided that its second funda- 
mental tensor b,; is of non-zero determinant, because, for such a Vj, 


Rigg = 5x) (5.2) 
(Eisenhart, 2, 197). Thus a classification of Riemannian spaces by 
means of the e-characteristic alone, without reference to the relation 
of the complex to the fundamental quadric, would treat as one class 
all hypersurfaces of a flat 5-dimensional space. This adds point to 
the remarks made in the last section about the significance of the 
g-characteristic. 

Stating the above result as a theorem, we therefore have: 

Every Riemannian V, which is a hypersurface of a flat V; is of 
e-characteristic [(111)(111)] at all non-singular points, provided that its 
second fundamental form is of non-zero determinant. 

If it is everywhere of g-characteristic [(111111)], that is, if R;,,; is 
a constant multiple of g;;,., then it is a space of constant curvature. 

When the determinant of a;; in (5.1) is zero, the Riemann complex 
is even more specialized. If the matrix [q,;] is of rank 3, the complex 
consists of all lines touching the cone a;;X‘X! = 0 and so is of 
e-characteristic [(222)] (Jessop, 3, 226; an example of this type of V, 
is given below). If [a;;] is of rank 2, the quadric a;; is a pair of planes, 
and the Riemann complex degenerates into the square of the special 
linear complex having the line of intersection of the two planes as 
directrix; that being so, the e-characteristic is now [(21111)] (Zindler, 
10, 1133; an example of this type of V, is also given below). The 
case when [q,;| is of rank 1 does not arise, because, if it were of 
rank 1, all components of the Riemann tensor would be zero. 

Although it is not proposed to find analytical conditions for all 
the special cases considered, it may be of interest to obtain by a 
geometrical argument a necessary condition for R;;,, to be of the 
form (5.1). The complex consists of tangents to the quadric a;,, 
which we assume for the moment to be non-degenerate. Hence the 
polar linear complex R;;,,p of any line p” is special, its directrix 
being the polar of p”/ with respect to the quadric. So R;;,,p™ satisfies 
the conditions for a line, whence 

fet R Rus = 0. (5.3) 
This is true for any line p”, and so must be a consequence of the 
identity 


ijmn 


— 
= 0 


2 
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satisfied by the coordinates of a line. Therefore 


where » is a scalar. Multiplying by «”"”¢ and summing for m, n, 


and hence, contracting r with p and q with s, 

(RR) = 6p, 
where (RR) = 


Substituting in (5.4), we deduce that, if the Riemann complex is of 
«-characteristic [(111)(111)], the Riemann tensor being of the form (5.1) 
with det \a;;| 0, then 

This result may also be established purely analytically. Reference 
to it was made in § 3 (3.32). 

Conversely, it can be shown that, if R;;,, satisfies (5.5), and if 
(RR) ~0, then the Riemann complex consists of tangents to a 
quadric. If, however, it satisfies (5.5) with (RR) = 0, then—ecf. 
Ruse, 5, 450 (4.4)—the complex degenerates into the set of lines 
meeting a conic, or dually, touching a cone (case [(222)]); or it 
degenerates still further into the set of lines meeting two different 
intersecting lines (case [(2211)]) or meeting two coincident lines (case 
[(21111)]). 

An example of the case of e-characteristic [(222)| is provided by 
the V, 

ds? = y, =z, = 2), 
for which = = Rag, = Rigis, 


the remaining effectively different components being all zero. In 
this case the equation of the Riemann complex at the point (2, y, z, t) 


(p+ = 0, 
which is the complex of lines touching the cone 
X?+Y?4+Z? = 0 


in S,. The Ricci quadric is degenerate and coincides with this cone, 
the only non-zero components of the Ricci tensor being 


Ry = —2 = Ry = Ry. 
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The g-characteristic of this space is easily seen to be [(111)(111)]. 
For the V, defined by 
ds? = 
all the components of the Riemann tensor are zero except 
= —P (= et.). 


When ¢ + 0, the e-characteristic is [(21111)]. The equation of the 
Riemann complex is _#(p2y? = 0, 
so the complex consists of the special linear complex p!* = 0 taken 
twice. The g-characteristic is [(11111)1]. It is perhaps rather sur- 
prising that a Riemannian V, can be as specialized as this without 
being actually flat. 

Reference was made in my previous paper (Ruse, 6, (6.13) et seq.) 
to the fact that any quadratic complex which is self-polar of the 
second kind with respect to the fundamental quadric is a harmonic 
complex. In general, if the Riemann tensor is of the form 


= Ay bj, — Og, (5.6) 


where a,;, b;; are symmetric, then the Riemann complex consists of 


ij? 
the lines which cut the quadrics a@;,,;, 5;; 


points, and is of the type known as harmonic. 


in harmonically conjugate 


The Riemann complex for every V, conformal to a flat space is 
harmonic, because, for such a Vj, 

Riga = — Gin (5.7) 
where d;; = $R,;—39;; R (Eisenhart, 2, 217-18). A harmonic com- 
plex is in general of e-characteristic [111111], but may have various 
other e-characteristics according to the relations between the quadrics 
a;;, b,; which determine it (Jessop, 3, 359). For example, the Rie- 
mann complex for the Schwarzschild space-time 

ds? = —(1—2m/r)-1 dr? —r? d6? —r* sin?0 +-(1—2m/r) dt? (5.8) 
is of e-characteristic [(11)(11)11] (Lamson, 4, 722), and is a special 
tj? bis 
of (5.6) have four common generators forming a skew quadrilateral. 
The singular surface (Ruse, 6 (5.15)) consists of two quadrics meeting 
in the same skew quadrilateral, their equations being 

= —p-*X2+ 2r? = 0, 


case of the harmonic complex in which the defining quadrics a 
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where » = 1—2m/r, while the fundamental quadric is 
943 = sin*9Z?+- pT? = 0. 
The quadrics a;;, b;; of (5.6) are in this case given by 


where w is a complex cube root of unity. Also 
= Bis, 

and the quadrics a;;, b;;, «;;, B;;, gi; therefore all belong to the same 
pencil and intersect in the same skew quadrilateral. Moreover, it is 
easy to see that the quadrics «;;, 8;; are polar to one another with 
respect to the fundamental quadric, and hence also that a,;, b,; are 
polar to one another. 

The g-characteristic of the Schwarzschild space-time (5.8) is 
[(1111)(11)]. In general, however, the g-characteristic of a harmonic 
Riemann complex, whether of the form (5.6) or of the more special- 
ized form (5.7), is [111111]. When the «- and g-characteristics of 
(5.7) are both [111111], there are no equalities among the e-roots 
Aj...) Ag or among the g-roots p44,..., 4g, but they are separable into 


= O= = 
(cf. Jessop, 3, 134). 
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ON LATTICE POINTS IN A CYLINDER 
By K. MAHLER $9(Manchester) 
[Received 18 May 1944] 


DENOTE by %,, Xp, #3 rectangular coordinates in three-dimensional 
space, and by K a convex body with the origin O = (0,0, 0) as its 
centre. A lattice A, . 


3 
Uz, (kh = 1,2,3; Ue, Ug = 0, F1, F2....), 


say of determinant 
d(A) = 


is called K-admissible if O is its only point which is an inner point 
of K. Let A(K) be the lower bound of d(A) extended over all 
K-admissible lattices. Then A(K) > 0, and there is at least one 
critical lattice, i.e., a K-admissible lattice A such that d(A) = A(K).* 
Minkowski’s theorem on convex bodies may be expressed as 


A(K) > 


where V is the volume of K. In general, the sign ‘ > ’ holds in this 
inequality, and so the problem arises of finding the exact value of 
A(K). Minkowski himself solved this problem for the cube, the 
octahedron, and the sphere. In this note I solve it for the cylinder 


K: —l<a<+l, 
by proving that7 A(K) = 43, (1) 


a result which surprisingly has escaped notice. 
That A(K) < $v3 is nearly trivial, because the following lattices of 
determinant }v3 are evidently K-admissible: 


(i) The lattices A, derived from the particular lattice 
= = = Us (a, arbitrary) 
by any rotation about the 2,-axis; 
* For the two-dimensional case of these rather obvious statements see my 
note, J. of London Math. Soc. 17 (1942), 130-3. 


+ The same proof shows that A(K) = 3v3, where K is the n-dimensional 
convex body 2{+23 < 1, |x5| < 1.,..., |a,| < 1(m > 2). 
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(ii) The lattices A, derived from the particular lattice 
= Uy, = (a, B arbitrary) 
by any rotation about the x;-axis.* 
It suffices therefore to show that 
A(K) > 4v3, (2) 
in order to prove the assertion (1). 
I use the following lemmas. 


Lemma 1. Let x be a plane convex polygon of area A with angles not 
greater than 120°, and let C,, C,,..., C, be non-overlapping circles of 


radius r contained in 7. Then 
A 


Lemma 2. Let n be a positive integer, and let W be the cube 


Let further Z,, Z,,..., Z, be non-overlapping circular cylinders of radius 
4 and height 1, all contained in W with their axes parallel to the x,-axis. 


Then 16 


t< 


Proof. Denote by x any number in the interval —n < x < n. The 
plane x, = x intersects the cylinders Z,, Z,,..., Z, in a certain point 


+n 
set J(x), say of area Q(x). Then the integral f Q(x) dx equals the 


total volume of the cylinders Z,, Z,,..., Z,, and so 
+n 
| Q(x) dx = at. (3) 
—n 


Now there are at most 2¢ different values of x for which the plane 
3 = x contains either the base or the top of one of these cylinders; 
let a be different from these exceptional values. Then J(x) consists 
of a finite number, say s, of circles of radius }; no two of these circles 
overlap, and all lie inside the square 
<2, <0, % = 2 

* The lattices A, and A, are the only critical lattices of K, as can be proved. 
One can further show that, if H is any convex body symmetrical in O which is 
contained in, but different from K, then A(H) < A(K). 

+ For a proof see the note ‘On the densest packing of circles’ by B. Segre 


and myself, American Math. Monthly, 51 (1944), 261-70. 
3695.17 c 
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of area A = 4n*. Hence, by Lemma 1, 
4n? 8 
= 
~(4)2v12 


Then 


and so, by (3), 


whence 


Proof of (2). Put 
so that the cylinder K consists of all points satisfying F'(x,,%,%3) < 1 
Denote by A any K-admissible lattice. Then at every point 
X = (x, v$, 78) of A place a cylinder 
Z(X): F(x,—2}, %.—2}, < 
of half the linear dimensions of K, and with its centre at R and axis 
parallel to the x,-axis. Since A is K-admissible and since K is convex, 
no two of these cylinders overlap.* 
Let x now be a large positive integer. Since every lattice parallel- 
epiped is of volume d(A), the cube 
< n—3, < n—3, < 
8n 
d(A) 
points X of A; at least as many cylinders Z(X) lie therefore in the cube 


contains 


+ O(n?) 


<0, <n, <n. 

Thus, by Lemma 2, 
8n8 
d(A) 


O(n?) < Sn’, 
V3 


whence 
d(A) > 4v3—o(1), ie. d(A) > A(K) > 43, 
as asserted. 


* Minkowski, Geometrie der Zahlen, 74. 


2 
+n 4 
int= [ Q@)de< | 
—n —n 
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ON AN INTEGRAL EQUATION ARISING IN THE 
THEORY OF DIFFRACTION 


By E. T. COPSON (Dundee) 
[Received 30 September 1944] 


1. Introduction 

THE problem of the diffraction of plane sound waves of small ampli- 
tude or plane electromagnetic waves by a perfectly reflecting half- 
plane was first solved rigorously by Sommerfeld* in 1896 when he 
introduced many-valued solutions of the equation of wave-motions 
having the geometrical shadow as branch membrane. Sommerfeld’s 
results were subsequently obtained by H. M. Macdonaldt who 
expanded the wave-function as a series of Bessel functions and by 
Lamb{ (for the case of normal incidence) who used parabolic 
coordinates. Although the subject already has a large literature,§ 
an interesting new contribution has recently been made by W. Mag- 
nus,|| who has shown that the problem can be reduced to the solution 
of a singular integral equation. 

Magnus arrives at his integral equation by the following ievsieal 
argument. If plane polarized electromagnetic waves with electric 
force E,e* in the direction of the axis of z are incident on the per- 
fectly conducting half-plane y = 0, x > 0, then, as Poincaré pointed 
out many years ago, an alternating current-sheet flows in the con- 
ductor in the direction of the axis of z and produces the scattered 
field #,e which is also polarized parallel to Oz. This field is 
given by 

0 
where the function f(€), so far unknown, gives the density and 
phase of the alternating current-sheet. But the total electric force 
(EZ, + £,)e™, being parallel to the perfectly conducting screen, must 
vanish on the screen. Hence f(€) must satisfy the integral equation 


[ dé = 
0 


* Math. Ann. 47 (1896), 317. + Electric Waves (Cambridge, 1902). 

t Proc. London Math. Soc. (2) 4 (1906), 190. 

§ See, for example, Baker and Copson, Huygens’ Principle (Oxford, 1939), 
chap. iv. \| Zeitschrift f. Phys. 117 (1941), 168. 


k 
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when x > 0, the function g(x) being the known expression to which 
—E, reduces on the screen. 

My interest in Magnus’s work was aroused because I had already 
obtained independently the same integral equation by a different 
method and in a somewhat different connexion, though at the time 
I was unable to solve the integral equation. My work was concerned 
with the diffraction of plane sound waves by a perfectly reflecting 
half-plane, or its analytical equivalent, the diffraction of plane 
polarized electromagnetic waves in which the magnetic vector is 
parallel to the edge of the screen. In this problem the boundary 
condition is the vanishing of the normal derivative of the velocity 
potential (or magnetic force) on the screen. By the argument given 
in §2 below I was led to Magnus’s integral equation with the im- 
portant difference—the screen is now the half-plane y = 0, x < 0. 

Magnus asserts that, on physical grounds, one would expect the 
current function f(€) to behave like a multiple of €-? when € is small, 
and he therefore assumes that 


=> Cn En ()- 


He then expands the function g(x) as a series of the form 


g(x) = An 


and obtains a set of linear equations for the infinite set of coeffi- 
cients* C,,. Having found the coefficients C,,, Magnus ultimately 
arrives at the Sommerfeld solution for normal incidence after long 
and complicated analysis. 

If this were the only method of solving the integral equation, the 
new solution of the Sommerfeld problem would have little merit— 
the introduction of the integral equation would in no way simplify 
the problem. If series of Bessel functions are to be used, it is simpler 
to use them from the beginning as Macdonald did. It tuins out, how- 
ever, that it is possible to give a direct solutior. of the integral equation 
by means of complex Fourier transforms and the theory of functions 
of a complex variable. This solution is given in §3 of the present 
paper, and the deduction of Sommerfeld’s formula follows in § 4. 


* The set of linear equations is one which had been solved rigorously by 
Titchmarsh, Math. Zeits. 25 (1926), 321, though Magnus appears to be 
unacquainted with this work. 
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2. The derivation of the integral equation 
Let us consider monochromatic sound waves of complex velocity- 
potential de”, where ¢ is independent of ¢ and w > 0. The time- 
factor e“ will be omitted hereafter. Then the velocity-potential ¢ 
satisfies the partial differential equation 
= 0, (2.1) 
where k = w/c and c is the velocity of sound. If plane waves with 
velocity-potential 
do — etkxcos a+ikysin « (0 <a< 


are incident on the rigid screen y = 0, x < 0, the total velocity- 
potential will be of the form 


= 
where ¢, represents the waves scattered by the obstacle. (See Fig. 1.) 
AN 


Direction of 
incident rays 


>x 


Fie. 1 


The function ¢, will be a function of x and y alone which satisfies 
the following conditions: 

(i) ¢, is a solution of (2.1) which is continuous and has continuous 
partial derivatives of the first and second orders in the (x, y)-plane 
supposed cut along the negative part of the x-axis; 

(ii) the normal derivative of ¢)+¢4, vanishes on y = 0, x < 0; 
(iii) at infinity, 4, satisfies Sommerfeld’s radiation condition.* 
These conditions determine ¢, uniquely. The radiation condition is 
essential; for otherwise it would be possible to add to the solution 

terms representing standing waves. 


* Sommerfeld, Jahresbericht der D.M.V. 21 (1912), 309. Cf. Baker and 
Copson, loc. cit., 25. 


Screen 
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In its original form the radiation condition is difficult to apply. 
A simpler form is to suppose that & is a complex number with small 
negative imaginary part so that 

—4n < argk <0 (2.2) 
and to require that ¢, vanish at infinity. Roughly speaking, this 
form of the radiation condition prevents ¢, from containing any 
terms proportional to e“*", so that ¢, represents only expanding waves 
at infinity. The condition (2.2) has a simple physical interpretation; 
for sound waves it means that the air has slight viscosity, whereas 
for electromagnetic waves, it asserts that the medium has slight con- 
ductivity. An additional advantage of this form of the condition is 
that, when k is complex, the Fourier-transform analysis is greatly 
simplified. 

Now a solution of (2.1) which is continuous and has continuous 
derivatives of the first and second orders in the half-space y > 0 and 
which satisfies there the radiation condition is given by the formula* 


7 on n=0 r 
where r2 = (w—€)?+ y?+ (2z—£)?, 
in the half-space y > 0. In y < 0, the sign of the right-hand side 
is reversed. By condition (ii), 


apy — —iksin f(z) (2.4) 


ey 


when y = 0, where f(2) is identically zero when x < 0 but is unknown 
when « > 0. Hence we have 


—tkr 
= F | | sin a #(€)}" 


=4hi {—ik sin « eos «+ f(¢)} dé, 


where R= 


* See, for example, Rayleigh, Sownd, 2 (1896), 107, equation (3). The 
formula is easily proved by the theory of Green’s function. In our problem, 
¢ happens to be independent of z. 


> 
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the upper or lower sign is taken according as y is positive or negative. 
It follows that 


= 1; f {—iksin dé. (2.5) 


If f(€) were identically zero, the problem would be the reflection 
of plane waves at a perfectly reflecting plane y = 0, and ¢ would be 


cos a+iky sina etkxcos a—ikysin « 


for y > 0, and zero for y< 0. In our problem f(é) is identically 
zero only for x < 0; hence (2.5) becomes 


etkxcos a+ikysin «| gikxcos a—tkysin a+ ER) dé (2.6) 
0 


for y > 0, and $= dé (2.7) 
0 


for y < 0. 

The velocity-potential ¢ is continuous across the aperture of the 
screen, that is, across the half-plane y = 0, x > 0. It follows from 
(2.6) and (2.7) that f(€) satisfies the integral equation 


[ dg = (2.8) 
0 


when a > 0. We have thus proved 


THEOREM 1. When plane sound waves with velocity-potential 
() << <m), are incident on the perfectly reflecting 
half-plane y = 0, x < 0, the total velocity-potential of incident and 
scattered waves is given by equations (2.6) and (2.7), where f(&) is the 
solution of the integral equation (2.8). 

In the problem considered by Magnus, plane polarized electro- 
magnetic waves, with electric force 


— pikxcosa+ikysin « 
_= etk. iky si 


parallel to Oz, are incident on the perfectly reflecting half-plane 
y = 0,x > 0. The integral equation for this problem can be obtained 
by a similar argument. For the electric force HZ, in the scattered 
waves is parallel to Oz and satisfies the following conditions: 

(i) Z, is a solution of (2.1), depending on x and y alone, which is 
continuous and has continuous derivatives of the first and second 
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orders in the (x, y)-plane supposed cut along the positive part of the 
axis of 2; 

(ii) H, = ony=0,z2> 0; 

(iii) #, satisfies Sommerfeld’s radiation condition at infinity. 
It follows that H, is an even function of y, 0#,/éy an odd function. 


Hence = as y>+0. 

But @£,/éy is continuous when y = 0, « < 0; therefore f(z) = 0 when 
x <0. Hence we have, by (2.3), 


—ikr 

0 -@ 0 
But since Z, = —H, when y = 0, x > 0, we again obtain the integral 
equation (2.8). 


3. Solution of the Integral Equation 
Magnus’s integral equation is of the form 


[ FQ HP dé = g(a), (3.1) 
0 


where g(x) = 2ie*reose when a > 0. Although this equation appears 
to be of the type = 
dé = (3.2) 


discussed by Titchmarsh,* it is actually quite different. In Titch- 
marsh’s work g(x) is given for all x, and f(x) has to be found for 
all 2; here g(x) is given for x > 0 and f(x) is identically zero for 
x < 0, and we have to find g(x) for x < 0 and f(x) for x > 0. 

In order to apply the theory of generalized Fourier integrals, it 
is necessary to make certain assumptions regarding the unknown 
function f(x). By (2.4), we have, when x > 0, 


= tksin f(x), (3.3) 
at li 
where ay 


Since —é@¢,/éy is the velocity parallel to Oy in the scattered waves, 
we might assume on physical grounds that (x) is continuous and 


* Titchmarsh, Theory of Fourier Integrals (Oxford, 1937), 314-15. 
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differentiable in 2 > 0, tends to zero as x > 00 by the radiation con- 
dition, and tends to infinity as x > 0 since the velocity is infinite at 
the edge of the screen. Magnus assumed that (2) becomes infinite 
like as x0. Actually it suffices to assume that belongs 
to L(0, 00) for any positive value of c. It follows that, if we write 
k=p—q (p,q > 9), (3.4) 
e~“*f(x) belongs to L(0,00) for > gcosa. 
If we define g(x) for negative values of x by (3.1), it follows from 
our assumption concerning ¥(x) that g(x) is continuous for x < 0 
‘and differentiable for x < 0. Moreover, since* 


H(z) ~ (=)te-«-, 


where |z| is large and —2z7 < argz < 7, we have 
HP (k|a—é|) = Of 


and so, when x = —z’, and 2’ is large and positive, 
0 


Therefore g(x) tends to zero as x > —oo, and also e“g(x) belongs to 
I(—oo,0) and to L?(—oo,0) for any positive value of c. This is 
what we should expect on physical grounds; for, comparing (3.1) and 
(2.7), we see that 

= 2% lim 


I next state certain results needed in the course of the solution of 
the integral equation. The complex Fourier transform of the kernel 
U(x) = H(k\a|) is denoted by 


L(w) = com j de — (-) Hi we de. 


Lemma 1. L(w) is an analytic function of the complex variable 
w = u+iv, regular in the strip —q < v < q, and has the form 


2 1 
where the branch of the square root reduces to k when w= 0. The 


* Watson, Theory of Bessel Functions (Cambridge, 1922), 198. 
+ Or even for c > —q, though this is unnecessary in what follows. 
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analytical continuation of L(w) is regular over the whole plane supposed 
cut radially outwards from w = k to infinity and from w = —k to infinity. 

To prove the formula for L(w), let us suppose in the first instance 
that |w| < |k|. Let us write ¢=argw, y= —argk so that 
0 < y< 4a. We may then rotate the path of integration through 
an angle —(47—y), to obtain 

wexp(— yi) 
V(37)L(w) = cos wa (ka) dx; 
0 
for, on the circular are z = Re~*, where 0 < 0 < 4n—y, we have 
|cos wz| < exp{R|w| . |sin(@—¢) |}, 
|H((ke)| = o[exp{—R|k| . |sin(0-+y) }], 

and so, by Jordan’s lemma, the integral along this arc tends to zero 
as Roo. If we now put ka = —it, we have 


0 


2 twt 
0 


J(k?—w?) 
by a known integral,* valid since 

| 9 
Im(*)| < 
< < 
The rest of Lemma 1 then follows by analytical continuation. 
The generalized Fourier integralst we shall use are functions of 

the complex variable w, defined by 


1 iwx 
F,(w) = | fleet de, 
0 


F_(w) = | da; 


and similarly for G,(w) and G_(w). In the present problem F_(w) 
is identically zero. 
* Watson, loc. cit. 388 (10). + Titchmarsh, loc. cit. 4. 


| at 
ik 
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Lemma 2. F,(w) is an analytic function of w, regular when v > q¢ cos «. 
It is of the form 
ksin « 
(277) w+k cos « 
where ‘¥(w) is regular in v > 0, bounded in v >c > 0. 
The formula for F',(w) follows from (3.3), with 


+¥(w), 


F,(w) = 


1 F 
¥w) = dex, 
0 


Since e-*%/(x) belongs, by hypothesis, to Z(0,00) for any positive 
value of c, this last integral converges uniformly with respect to w 
inv >c. Hence Vw) is regular in v > 0. Moreover, if v >, 


l 
[F(w)| < | dex. 
v(2z) 


3. 


G(w) = when v > qcosa. 


w+k cos « 


This is evident, since 
g(x) = DMetkrcosa 

when a > 0. 

Lemma 4. G_(w) is an analytic function of w, regular in v < 0, 
and tends to zero as u-> +00 for any fixed negative value of v. 

Since e“g(x) belongs to L(—oo, 0) for any fixed positive value of c, 
the integral defining G_(w) converges uniformly with respect to w in 
v < —c. Hence G_(w) is regular in v < 0. Moreover, 


0 
G_(u—ic) = com J da, 


which tends to zero as u > +00, by the Riemann-Lebesgue theorem.* 
5. > qceosa, b > 0, 


o+ai 
x) dw G_ dw 
in the mean-square sense. 
* Titchmarsh, loc. cit. 11. 
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The function h(x), which is equal to e-**g(x) when x > 0, zero 
when x < 0, belongs to L*(0,00). Its Fourier transform is 


H(u) = em J dx = G,(u+ia). 


| du 


Hence h(x) = —— 

in the mean-square sense. Similarly, if k(x) is equal to e°%g(x) for 
x < 0, zero for x > 0, 


k(x) = | du, 


again in the mean-square sense. The result follows by addition. 


Lemma 6. If gcosa<a<q, 
f(QHP(k\a—E|) dé = (w) L(w)e-#” dw 
0 —ao-+ia 
in the mean-square sense. 
The function h(x) = e-**H(k|x|) belongs to L*(—oo,0o) if 
—q <a<gq. Its Fourier transform is 


] . 
H —axrpiux ia). 
(u) = | H® du = L(u+ia) 


The function k(~) = e-f(x), where f(x) = 0 when x < 0, belongs 
to L(—o, 00) if a > qceosa. Its Fourier transform is 


K(u) = dx = F,.(u+ia). 


The result now follows from a known theorem.* 


Lemma 7. Let H(w) be regular in the strip a,Svu<a,. Let 
H(u+iv) be L(—o0,00) or L*(—0co,00) and tend to zero as u—> +00 
for v in the above strip. Let K(w) have similar properties in b, <v < bg, 
where b, <a,. Let 


o+ia o+ib 
[ dw = dw (3.5) 
—o+ia —2-+ib 


* Titchmarsh, loc. cit. 90, Theorem 65. 
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for all x, where a, << < dy, <b < b,. Then H(w) and K(w) are 
regular and equal in the strip b, < v < a, and tend to zero as u—> +00 
uniformly in any interior strip. 

This is a known result.* 

We now apply these lemmas to the solution of the integral equation 


dé = g(x). 
0 


By Lemmas 5 and 6 we have, if gcosa <a <g,b>0, 
o+ia co —ib 
dw = G_(w)e-*” dw 

—o-+ia 
in the mean-square sense. 

In the strip gcosa+e < v < q—e (ec > 0), the function 

H(w) = 
is regular, and H(u+iv) tends to zero as u > +00 for v in the above 
strip; this follows from Lemmas 1, 2, 3 and the fact that Y(w) is 
bounded. Moreover, H(u+iv) belongs to L?(—0o, 00) for any field v 
in the strip. 

In the strip —b’ << v < —e (ce, b’ > 0), G_(w) is regular and 
G_(u+iv) tends to zero as u > +00 for v in the strip. Lastly, since 
e%g(x) belongs to L?(—oo, 0) for any positive c, G_(u+-iv) belongs to 
L?(—0o, 0) for any fixed v in the strip. 

The conditions of Lemma 7 are thus all satisfied. It follows that 
H(w) and G_(w) can both be continued analytically throughout the 
strip —b<v<qg and are equal and regular there; moreover, 
H(u+iv) and G_(u+iv) tend to zero uniformly as wu — +00 in any 
interior strip. 

But G_(w) is regular in v < —6 for any positive value of 6. Hence 
we can continue H(w) all over the half-plane v < g, and so 

= @_(w) 
when v < q. 
By Lemmas | and 3 we now have 


Lemma 2 shows that F,(w) is regular in v > 0 apart from a simple 
pole at w = —kcosa. By (3.6), F,(w) can therefore be continued 


* Titchmarsh, loc. cit. 255, Theorem 141. 
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analytically all over the w-plane, and the only singularities of the 
function so defined are the above simple pole and a branch-point at 
w = k in the lower half-plane. We may therefore write 


(3.7) 


w 


where ¢(w) is an integral function. 
Substituting this value of F,(w) in (3.6), we obtain 


= — (=) a)G(w)|. (3.8) 


Since G_(w) is regular at w = —k cosa, it follows that 
(k—k cos «) 


+(w+k cos «)x(w), 


¢(—kcosa) = — 
V(k—k cos «) 


where x(w) is also an integral function. Hence we have, by (3.7) and 
(3.8), 


and so 


Kk—k cos x) ./(k—w) 


aw) = 1 4 2x(w) (3.10) 


J(k+w) Jw+keosa 
It remains to determine y(w). 

Now F,,(w) is bounded in v > c, where gcosa < c < q; hence, by 
(3.9), x(w) is O(w-*) as |w| >co inv >c. Again G_(w) is bounded 
in v < —c, and so, by (3.10), y(w) is O(w) as |w| +00 in v < —c. 
Lastly G_(w) tends to zero as w-> +o uniformly in —e < v <¢; 
hence x(w) is o(wt) as |w|-—> oo in this strip. Hence the integral 
function x(w) is O(w?) as |w|-—>0o. By the extension of Liouville’s 
theorem* x(w) is a polynomial of degree < }, and so is a constant; 

and, as x(w) is O(w-) in v > ¢, this constant is zero. Hence we have 
(k—k cos «) ./(k—w) 
(27)(w-+-k cos x) 


We have thus proved 
THEOREM 2. Let k = p—igq, where p>0,q> 0. Let f(x) be of 
the form f(x) = iksin a «4 
* Titchmarsh, The Theory of Functions (Oxford, 1932), 85. 
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where e-%p(x) belongs to L(0,00) for any positive value of c, and 
0<a<a. Then the integral equation 


[ HP dg = 0) 
0 


has the solution whose complex Fourier transform is 


Fle) = | dz = — 008 a) (3.11) 


when v > qeosa. The explicit form of the solution is 


V(k—k cos «) (k—w) 
f(x) = ak sin « 5 e dw. (3.12) 


Although it is not necessary to find f(x) explicitly in the solution 
of the diffraction problem, it is nevertheless of some interest to give 
the actual solution of the integral equation (3.1). We have 

| F,(w)e*™” dw 


1 


f(x) = 


iat 
\(k—k cos x) ./(k—w) dp 
2a w+k cos « 


where gcosa <a<q. Deforming the path of integration into the 
real axis, we obtain the required formula. 

If we make the imaginary part of k = »p—iq tend to zero, the 
equation (3.12) will still hold, provided that we indent the real axis, 
the indentation being downwards at w = —kcosa and upwards at 
w = &, 


4. The solution of Sommerfeld’s problem 

Instead of solving the diffraction problem by using equations (2.6) 
and (2.7) to get ¢(z,y) from f(x), it is much simpler to make use of 
the transform F,(w) in the following way. 

When y < 0, we have, by (2.7), 


$(x,y) = —¥ | fle)m(a—€) dé 
0 


+ia 
F,(w)M(w)e-*” dw, 
+ 


—o-+ia 
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where m(x) = AP {k 
M(w) = : m/(x)e*” dx 
al(2z7) 

and qceosa << a< q. 


The function M(w) is regular in the strip —q < v < q; its value is* 


M(w) = 
tr} 


where the branch of the square root reduces to k when w = 0. The 
analytical continuation of M(w) is regular in the whole w-plane sup- 
posed cut radially outwards from & to infinity and from —k to 
infinity. Thus ¢(2,y) is determined for y < 0; and similarly for 
y > 0. The resulting solution of Sommerfeld’s diffraction problem 
is contained in the following theorem. 


THEOREM 3. When plane sound waves with velocity-potential 
etkxcosatikysina (() << aq < m) are incident on the perfectly reflecting 
half-plane y = 0, x < 0, the total velocity-potential of incident and 
scattered waves is given by 


d(x, y) — eikxrcos a+ikysin gikxcos a—ikysin «| 


o+ia 
+ V(k—k cos a) dw 
Qri V(k+w) w+keosa 
when y > 0, and by 
o+ia 
d(x y) = — V(k—k cos a) dw 
Qari J(k+w) w+kcosa 


when y < 0, where qcoosa<a<q. 

To identify this form of the solution with that given by Sommer- 
feld, let us consider the case when (x,y) lies inside the geometrical 
shadow, so that 


x = —pcos8, y = —psin8g, 


* This result is the particular case of Watson’s equation—loc. cit. 416, 
equation (2)—with = 
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where p = +./(2?+y?), 0< @< a. By Cauchy’s theorem the path 
of integration can be deformed into the hyperbolic path 


w= —kecos(0+ir) 


without crossing the pole w = —kcosa or the radial cut from —k 
to infinity. (See Fig. 2.) This gives 


= {1—cos(6-+-i7)} cos(@+-i7)—cos « 


—ikpcosht 


—ikpcosht 1 1 


1 1 


] 
= cosh t—cos(a—@) cosh r—cos(a+-@) 


0 


— — | e-tkpcosht cos dy cosh $7 cos 38 cosh dr 
7 cosh7+cosy cosh7-+-cos8 
0 


where 8 
3695.17 D 


ig 
ag 
x 
Fic. 2 
| 
1 
1 
0 


34 ON AN INTEGRAL EQUATION 


This can be transformed into* 
4(2kp)sin 4(0—«) 


ikpcos(@—a) 
y) = + 
—*(2kp) sin 3(0+«) 
+ e-tkpcos(0+a) al 


which is Sommerfeld’s result. The other cases can be handled in 
a similar manner. 


* See, for example, Baker and Copson, loc. cit. 139-40. 
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PROJECTIVE SPACE 
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Let @, @’ be two curves in a projective n-dimensional space S,,, and 
suppose that @ and @’ touch each other at a (possibly singular) 
point O, at which they have the same osculating spaces of dimensions 
1, 2,..,r(1<r< n—1). Whenr = n—1, @ and @’ have (at O) n—1 
projective tac-invariants 4, J4,,..., %,_,, of which some simple geo- 
metric interpretations are known.* This has been extended to the 
case in which 1 < r < n—1 by B. Su,t who has proved that @ and 
@' have then (at O) r—1 projective tac-invariants 

(1) 
This extension can be deduced almost immediately from the previous 
result, by considering the projections @,, @; of @, @’ upon their 
common osculating S, at O from an arbitrary S,_,_, skew to §.. For 
@, and @; have at O the same osculating spaces of dimensions 
1, 2,..., r—1, and so the previous result can be applied. The r—1 
projective tac-invariants, (1) say, thus arising from @,, @; at O, are 
easily seen to be independent of the centre S,,_,_, of projection, so 
that they also are projective tac-invariants of @ and @’ at O. 

The last result gives nothing if @ and @’ touch at O but have 
distinct osculating planes at O, i.e. when 1 = r < n—1. In this case, 
however, C. C. Hsiung has very recently obtained by direct calcula- 
tion a number of projective tac-invariants, on the hypothesis that 
both @ and @’ have at O an ordinary simple point. Moreover, he 
has given an interesting, but rather intricate, geometric interpreta- 
tion of these invariants, involving certain irrationalities. These 
results are now extended and completed by the following 


TuroreM. Let @, @' be two curves in S, = S), (n > 3) having an 
ordinary simple point O in common. Denote by S,,, Sj, their respective 
k-dimensional osculating spaces at O (k = 1, 2,..., n—1), and suppose 
that .., (2) 

* Cf. B. Segre (2) also for further references on the subject. 


+ Cf. B. Su (3). I know this paper only through the review appearing in 
Math. Reviews, which contains no indication about the proof. 
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where r is any fixed integer satisfying 1< r < n—2, but that no other 
particularization of these spaces occurs. Hence, if r+1<i< n, the 
intersection of S;_, and S},,,_; is S,, and not a space of higher dimension, 
so that the join of S,_, and S},,,_; is a prime, which we indicate by p;. 

Let us fix in S,, an arbitrary prime 7 not containing O, and further 
consider a variable prime x not containing O and tending to a position 
Xo through O, but not through S,; hence x meets the line S, at a point O, 
tending to O. The prime x has then exactly one point in common with 
each of the curves ©, @' in the neighbourhood of O. We denote these 
two points by Q, Q’, and the intersections of their join with the primes 


a, p; by P, R; respectively (¢ = r+-1, r+-2,..., ”). 
When x tends to xo, the expression 
has a limit, which we denote by 4,; (i Aj; =7r+1, r+2,..., 


this limit is independent of 7, xo, as well as of the way in which x tends 
to xo. Each of the numbers 4,; thus defined is consequently a projective 
tac-invariant of € and @' at O. 

The tac-invariants J4,; satisfy the identities 


and, when n > r+3, 

SI xX x 1, (5) 
where i, j, h are any three different numbers out of r+1, r+2,..., n. 
Hence, in any case, each of the tac-invariants 4,; can be expressed as 
(6) 
The n—r—1 tac-invariants (6) are independent; moreover, when r > 1, 
they are also independent of the r—1 tac-invariants (1) previously 
considered. 

When r = 1, the tac-invariants 4; coincide with those otherwise 
defined by C. C. Hsiung (1), who has dealt only with this particular 
case. It may be added that the geometric interpretation of the tac- 
invariants 4,; given above, is remarkably simpler than that obtained 
by this author for r = 1. 

Proof. Let us introduce in S, any system of non-homogeneous 
projective coordinates 2, Xp,..., x, with origin at O, having the axes 
of x, %,..., X, as independent lines lying in the r spaces (2) respec- 
tively, and the axes of 2,,,, a8 independent lines not lying 
in S,, but lying in the (r-+1)-dimensional spaces of intersection of 


Ij, = 1, (4) 
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and and S}_,,..., 8, and S},, respectively. Then @ and 


@’ can be represented parametrically by the equations 
2, = t, % = {a,t%...., 7, = {4,0}, 244 = = {at} (7) 
and 
= 0, = = {40%}, 245 = 
(8) 
respectively, where ds,..., @,, @,, are 2n—2 non-zero in- 
dependent constants, and the brackets { } denote that we have 
written only the terms of lower infinitesimal order in ¢ or t’. When 
ry > 1, the r—1 tac-invariants (1) of @ and @’ at O are 
J, J, Fy = a,/a,,. (9) 
It is immediately seen that the prime p;, previously defined for 
¢=r-+l1, r+2,..., n, has the equation 


x; => 0. 
Moreover, we can take the equation of z in the form 
C1 2%, = 1, (10) 
and the equation of x in the form 
+ =€ (€ #9), (11) 


since neither 7 nor x contains the origin O, and x meets the 2,-axis 
at a point O,(e, 0,...,0) lying in the neighbourhood of O, so that all 
the coordinates of this point must be finite. When x varies, the 
coefficients €, y,..., Y, in (11) also vary; and, when x tends to xo, 
these coefficients have finite limits, the first of which is zero. The 
points Q and Q’, i.e. the intersections of y with @ and @’ respectively 
which lie in the neighbourhood of O, are given by (7) and (8), where 
the parameter ¢ or ¢’ has to be chosen in the neighbourhood of 0 in 
such a way as to satisfy (11); hence 
lim = 1, lim ¢/t’ = 1. (12) 
We notice now that, when A varies, the prime 
Lz = 0 (r+1<ign) (13) 
describes a pencil, containing p; and 7 for A = 0 and A = 00 respec- 
tively. The primes of this pencil passing through Q and Q’ arise 
from certain values of A, say A = p» and A = yp’ respectively; from 
(7), (8), (13) we see that 


p= {a; p’ {a; 


| 
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where the brackets { } have the same significance as before. It follows 
that 
(QQ'R; P) = (up’0 00) = p/p’ = {a,th 

since the cross-ratio of the points Q, Q’, R;, P equals the cross-ratio 
of the primes of the pencil (13) containing them, and so (3) gives 

The right-hand side of this equation remains unaltered if we multiply 
the expressions in square brackets by and 
en ir—2i+1 — respectively. Then, from (12), we see that, 
when x tends to xo, 4;; has in fact the limit 
= "47-2841, (14) 


which does not depend on the coefficients of the equations (10), (11). 
The projective invariance of 4,; is now an immediate consequence 
of the projective invariance of the cross-ratio. Moreover, (14) shows 


that the tac-invariants 4; satisfy the identities (4), (5), and that 


the n—r—1 tac-invariants (6) are independent, because they have 
a non-zero Jacobian determinant with respect to the n—r—1 argu- 
ments @,- Finally, if r > 1, the tac-invariants (6) are 
independent of the tac-invariants (1), since the expressions (9) for 


the latter involve none of these n—r—1 arguments. 
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1. Introduction 

Ir is known that, if two plane curves touch at a point O, having 
contacts of order k—1 with their common tangent at O, and if non- 
homogeneous projective coordinates (x, y) are so chosen that the 
power-series expansions representing these two curves in the neigh- 
bourhood of O are y = az*+... and y = ba*-+-..., then the ratio a/b 
is a projective invariant of the curves, which was first found by 
H. J. 8. Smitht and R. Mehmkef in the case k = 2. Next C. Segre 
gave§ the invariant in the general case, with a simple geometrical 
characterization, and extended|| the case k = 2 to two curves in 
ordinary space touching at a general point with common osculating 
plane. Then B. Segre] further extended this case to two curves in 
a projective space S,, of n dimensions having at a common point the 
same osculating linear spaces of dimensions 1,...,»—1. Recently, as 
a supplement to B. Segre’s investigation, B. Su** has studied two 
curves in space S, having at a common point the same osculating 
linear spaces of dimensions I,..., k (1 << k < n—1), and characterized 
some projective invariants of the curves, using two theorems due to 
S. 8. Chern.tf 


+ H. J. S. Smith, ‘On the focal properties of homographic figures’: Proc. 
London Math. Soc. (1), 2 (1869), 196-248. 
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It is the purpose of the present paper to supplement the investiga- 
tions of B. Segre and Su, by studying two curves in S, having only 
one point and the tangent at it in common. This investigation may 
be regarded as a generalization of another one by the author.t 


2. Derivation of invariants 

Let C,, C, be two curves in a projective space S,, touching at a 
general point O, and suppose that the osculating plane of either of 
them at O is not contained in the osculating hyperplane of the other 
at O. Further, let 2,,..., x, represent non-homogeneous coordinates 
of a point in S,. If we choose O as origin, the common tangent of 
C,, C, at O as x-axis, the axes of x, and z, in the osculating planes 
of C, and C, at O respectively, and the x,-axis (¢ = 3,..., n—1) in 
the plane of intersection of the osculating linear space of ¢ dimensions 
of C, at O with the osculating linear space of n—7+-2 dimensions of 
C, at O, then, in the neighbourhood of O, the equations of the two 
curves can be written in the form 


Ci: = = 2,..., 2), (1) 

Cy: = **7+... (8 = 2,..., 0). (2) 

To find projective invariants of contact of C,, C, at O, we consider 

the most general projective transformation of coordinates which 

leaves the point O and the common tangent invariant, and changes 

the other axes in such a way that the new axes have the same pro- 

perties as the old. This transformation is expressed in terms of the 
non-homogeneous coordinates by the equations 


n n \ 
k=1 K=1 
* n 
= [(i+2 (2 = 2,..., 2) 


The effect of this transformation on (1), (2) is to produce two other 
systems of equations of the same form whose coefficients, indicated 


(3) 


by asterisks, are given by the formulae 
Vig = | 


(4) 


* — 9 
OF = ‘* bi n-i+2 = 2,..., J 


+ C. C. Hsiung, ‘Some projective invariants of certain pairs of space 
curves’, to be published in America. 
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On eliminating the a’s from the equations (4), we see immediately 
that the expressions 


n—2j+2 a n—2i+2 
I; = ‘| (i $j; i, = 2,...,) (5) 
Di Bs n—j+ 


are projective invariants associated with the point O of contact of C,, C,. 
In particular, if 7 = n—i+2 and if the space S,, is of odd dimen- 
sions, that is, if m = 2m+-1, m being a positive integer, then from 
(5) we easily see that the expressions 
Din—i+2 
are invariants. If, on the contrary, the space S,, is of even dimen- 
sions, that is, if m = 2(m-+-1), then besides the invariants (6) there 
is the additional invariant 


(7) 


Tn, +2,m+2 = b 
m+2,m+2 


3. Geometrical characterizations of the invariants J;,,_;,. 

Let the ‘point at infinity’ on the 2;-axis (¢ = 1,..., n) be denoted 
by O;. Then, from §1, we first obtain the following geometrical 
characterization: 

The invariant Jn.2m+2, associated with the point O of contact of the 
curves C,, C, immersed in a space Sy,,,1, is Segre’s invariant of contact 
at O of the plane curves obtained by projecting C,, C, from the space 
On to the plane OO, O,, 

In the second place, we shall characterize geometrically a general 
invariant J; ,_;,. at O of the curves C,, C, immersed in a general space 
S,,. For this purpose we may assume, without loss of generality, that 
2<%< n—i-+2, and consider the principal planey at O of the two 
projections Tj, [, of C,, C, from the centre 

O,...0;-4 On-i+s---On 
on to the space OO, O; O,,_;,.. It is immediately seen that the equa- 
tions of the space OO, O; O,_;,. are 

Wy = = = = = = = O, (8) 
and that the equations of I, [’, are (8) and 

Ty: = = On Mit. (10) 

+ G. H. Halphen, ‘Sur les invariants différentiels des courbes gauches’: 

J. de U’Ecole Pol. 28 (1880), 25. 
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To obtain the principal plane we project [, [, on to the plane OO, O; 
from a point V arbitrarily chosen in the space OO, 0;0,_;,5. If 


are the coordinates of V, the two projections are easily found to be 
given by equations (8), 2,_;,. = 0 and 


(12) 


n—-i+2,i = 

respectively. The cones projecting |, [, from V have contact of 

order at least i along OV if, and only if, the invariant of contact 

of the curves (11), (12) at O equals 1; then, and only then, the centre 

V of projection lies on the principal plane. Hence this plane has the 

equations (8) and 

= O. (13) 

We now consider the system of cones of order n—i-+-1 and vertex 

O in the space 00, 0;0,_;,. such that the polar surfaces of orders 

n+2i+1, n—2i+3,..., n—i+1 of any point in the plane 

OO; O,,-;+9. With respect to any cone (the last of these polars being 

the cone itself ), all pass through the line OO,. The equations of this 
system of cones can easily be written as (8) and 


J, k=0,1,...,t- 
j+k=i-1 j+k=n—-it+1 


th = 0. (14 


If the cones of this system have a contact of the i(n—7i-+1)th order 
with T, at O, then their equation (14) reduces to 


i-1 yn—i+l ak 
vy An —i+2,n-i+2 Vi + jk i n—i+2 0, 


n—it+2 
(15) 
where the A’s are arbitrary constants. Moreover, we determine other 
cones of this system by the conditions that they have contact of 
order i(n—i+1) with [T, at O. The equations of these cones are 
likewise found to be (8) and 


~n-2i+2 ~i-1 i—1 n—i+1 ak 
(16) 


where the B’s are arbitrary constants. The cones (15), (16) intersect 
in n lines residuaily to OO,. The planes projecting these m lines from 
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OO, are represented by the equation obtained on eliminating 2, 
from (15) and (16); they may reduce to a single plane to be counted 
n times, which then coincides with one or other of the n planes 
represented by the equations (8) and 
n— OR i% -i+2 = 0. (17) 
Thus we obtain the following geometrical characterization of the 
invariant J; »_;+9: 

The cross-ratio of the planes OO, O,,_;,2, OO, O;, (13), and any one 
of the n planes (17), is equal to one of the values of —'NJiz} ;.>. 

This characterization is simply an extension of the one already 
given for the invariant of two curves in ordinary space. 


4. A geometrical characterization of a general invariant /,,; 

Finally, I shall give a geometrical characterization of a general 
invariant J;; of the curves C,, C, of S,. Without loss of generality, 
it may be assumed that 2<i<j. On projecting the curves C,, 
C, from the space 0O,...0;_,;0;,;...0;-10;.;...0, on to the space 
O00, 0; 0;, we obtain two curves K,, K,, given by the equations 


and 

Ky: % = gag = (20) 


Firstly, as in the previous section, we consider in the space OO, O; O,; 
a system of cones of order j—1 and vertex O, having contact of 
order i(j—1) with the curve K, at O and satisfying the conditions 
that the polar surfaces of orders 1,..., j—i—1, j7—i+1,..., j—1 of any 
point in the plane 0O;0; with respect to any of them, all pass 
through the line OO,. The equations of this system are easily found 
to be (18) and 


where the £’s are arbitrary constants. Similarly, we have in the 
space OO, O; O; another system of cones of order n—i+ 1 and vertex 
O, having contact of order (n—i+1)(n—j+-2) with the curve K, at O 


+ C. C. Hsiung, loc. cit. 


m=1,...,3—1 
l+m=j-1 
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and satisfying the conditions that the polar surfaces of orders 
i... 
of any point in the plane OO;0; with respect to any of them, all 
pass through the common tangent OO,. The equations of this 
system are (18) and 


“jie 


j—t—1, n—1+1 


j—t 
in—i+2% 


l+m=n—i+1 
where the F’s are arbitrary constants. If we further impose on the 
cones (21), (22) the conditions that the planes through OO, in which 
their intersection lies reduce to a single plane to be counted n times, 
then these planes coincide with one or other of the » planes satisfying 
the equations sf and 


Next we issues in the space OO, O; 0, any algebraic surface of 


order n(j—7)+n—i-+1, satisfying the follow ing conditions: 


(i) the polar surfaces of orders 
1,..., (n+1)(j—t)—1, (nN 


of any point in the plane OO, 0; with respect to the algebraic surface 


all pass through the point 0,; 
(ii) the polar surfaces of orders 
1,..., n(j—t)—1, n(j—t)+1,..., 
of any point in the plane O, O; O, all pass through the point O; 

(iii) the curve of intersection of the surface with the plane O, O; 0; 
consists of n(j—i)+n—i+1 lines of which (n+ 1)(j—t) coincide 
with O;0,, and the other n—j+1 coincide with O, O;. 

It is easily seen that the equations of any of these algebraic sur- 
faces may be written as (18) and 

ain > tay = @, (24) 
We uniquely determine a surface in the system (24) by the con- 
dition of having a contact of the j(n—i+1)th order with K, at O. 
From (19) it follows at once that the equation (24) reduces now to 

Similarly, we uniquely determine another algebraic surface of 
order n(j—i)+j—1 in the space OO, O; O,, satisfying the following 
conditions: 
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(i) the polar surfaces of orders 
of any point in the plane 00,0; with respect to the algebraic 


surface all pass through the point 0O,; 
(ii) the polar surfaces of orders 
1,..., n(j—t)—1, n(j—1)+1.,..., n(j—t)+j—1 
of any point in the plane O, O; 0; all pass through the point O; 

(iii) the curve of intersection of the surface with the plane 0, 0; 0; 
consists of n(j7—i)+j—1 lines of which (n+-1)(j—7) coincide with 
O;0;, and the other i—1 coincide with O, O;; 

(iv) the surface has contact of order (j—1)(n—i+2) with K, at O. 

The equations of this surface are found to be (18) and 


The planes through OO, in which the curve of intersection of the 
two algebraic surfaces (25), (26) lies, are given by the equations 
(18) and 


n—it+1pi-1 


Thus we arrive at the following geometrical characterization of a 


general invariant J;;: 

The cross-ratio of the planes OO, O;, OO, O;, any one of the n planes 
(27), and any one of the n planes (23), is equal to one of the values 
of 
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NOTE ON THE DISTRIBUTION OF THE INTERVALS 
BETWEEN PRIME NUMBERS 


By LORD CHERWELL (Ozford) 
[Received 29 April 1945] 


In this note an attempt is made to set out some of the results con- 
cerning the frequency of prime-pairs, -triplets, and so forth which can 
be derived by probability methods from the assumption that the 
distribution of the prime numbers may be treated as ‘random’, This 
approach clearly precludes any attempt at rigorous proof in the 
mathematical sense. Nevertheless, it seems interesting that results 
obtained by such simple means fit fairly closely the facts found by 
our (necessarily limited) enumeration from tables, and that in two 
cases at least formulae emerge which are identical, as has been 
pointed out to me, with those derived by more elaborate methods.* 

If P is an odd prime number P-+-n can only be prime if n is even. 
In what follows P and p always denote odd primes (p ~ P) and n 
always an even number. 

If we assume a random distribution of primes, it is easy to show 
that the probability of P+-n being a prime is enhanced by a factor 
(p—1)/(p—2) for each different odd prime p which divides n. Let us 
compare for instance the probability of P+2p not being divisible 
by p with the probability of P-+-2 not being divisible by p. If we 
know nothing about P, the probabilities are clearly equal, namely, 
each (p—1)/p. But, if we know that P is a prime, then neither P nor 
P-+-2p can be divisible by p.t On the other hand, since one of the 
p—1 intermediate odd numbers must be divisible by p, the chances 
of its not being P+ 2 are (p—2)/(p—1). As the chance of P+-2p not 


* T have to thank Professor E. Maitland Wright not only for criticizing this 
note and suggesting many improvements but also for drawing my attention 
to previous work in this field, viz. Hardy and Littlewood, Acta Math. 44 (1923), 
1-70, and Staeckel, Sitz. der Heidelberger Akad. der Wiss., Abt. 1916, No. 10, 
1-45; 1917, No. 15, 1-52; 1918, No. 2, 1-47, and No. 14, 1-67, which also 
contain references to earlier work. A note on the connexion of this work with 
mine is at the end of this paper. 

+ If p were equal to P, this would of course not be true. As in our con- 
siderations P may always be taken as greater than }n, it cannot be a divisor 
of n. In any event the correction would be insignificant unless P were very 
small. 
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being divisible by p is 1, the probability, so far as the factor p is 
concerned, of P+-2p being prime, if P is prime, is 
1 
= 14+—— 
p-2 
.times greater than the probability of P+-2 being prime. Since this 
reasoning holds for any value or multiple of p, the formula is quite 
general, and we may write for the probability of P+-n being a prime 
if P is one, W, = W. [] (p—1/(p—2), the product being taken over 


p>3 


all different odd prime divisors of n and W, being the probability that 
P-+-2 is a prime. It is of course clear from this formula that W, has 
a minimum value W, when n is any power of 2. 

In order to test this formula we may count how often the various 
intervals from 2 to, say, %) occur in a given region in the neighbour- 
hood of N. According to the prime number theorem there will be 
about AN/log N primes in the region AN. If, is small compared 
with AN but large compared with log NV, there will be between each 
of these primes P and P-+-n, about n,/log N primes. In the region 
AN therefore we shall count a total of something like n)» AN/(log NV’)? 
intervals. If the probability of finding an interval n is 


the total number of intervals counted must be equal to 


AN 1 
W, 1 
n=2 pl\n,p>3 
As there are 4n, even numbers between 0 and mp», this sum will 


have $n, terms, so that, if we denote the mean value of 


1 
p—2 


for all even intervals from 0 to vo by s,,, we have the equation 


No Sng AN 2 AN 
8, log N’ 
2 
log 


no 


from which it follows that W, = 
pin, p23 


i.e. 
l 
. 
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The mean value s,,, of 


pin pin 
between 0 and n, may of course easily be found by summing all the 
}n, terms and dividing by 3n). We have 


n n 2 

—°! terms 1 contributing to the average —° xlx—=l1;: 

2 n 
= 0 


] 
FA terms 5 contributing to the average 
2p —2 
No 1 


p(p—2)’ 


1 
Een terms So Contributing to the average a further 
Pe (P1—2)(P2—2) 
My 1 2 1 


\ 


and so on. 
Hence, ignoring the square brackets, it is plain that s,, converges 


to a limit for large values of no. 


= 1-5147...* 


each prime factor, of course, being counted only once so that p runs 
from 3 to 00, p, is always greater than p,, etc. 
* It may be shown that the error is less than (log 4% 9)4/m) where A is some 
suitable constant. For 


1 2 1 
ln > + — 


Di, 


2 
2 2(1 
< — exp(A loglog = 2(log 
No 


The error clearly tends to zero as n)>—> 00. For the introduction of loglog 4mq, 
see Hardy and Wright, Introduction to the Theory of Numbers, Chapter 22, 
Theorem 422 or 429. 
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In the limit therefore the number of prime pairs of the form P, 
P-+2,—and for that matter P, P+ 4, or P, P+2”,—in an interval 
AN will be 
2 
S,log 
if Aw is the number of primes in the interval. If we put 


_ AN 

jog N 

_  2AN_ _ 1-320...AN 
S,(logN)2? (log V)? 

The accuracy with which this formula—which is identical with that 
derived by Hardy and Littlewood and by Staeckel—represents the 
facts is altogether surprising. Indeed, as is shown in Table 1, it is more 
exact than the mere probability approach gives us any right to expect. 


An 


An 


this becomes F, 


TABLE | 
| P,+F,/2 
AN.10-* | An(obs.) | W, (cale.) | F, (cale.) | (F, obs.) | F, (obs.) (obs.) 
220 to 255 | 2,807 0-1065 | 2995 | 295 | 299 297 
90 to 110 1,735 | 0-1145 | 199 | 198 | 206 202 
10to 30 | 2,004 | 0-1333 | 267 | 262 | 267 264-5 
5to 10 562 | 0-1480 | 832 | 79 82 80-5 
2to 5 | 364 | 01619 59 62 
O5to 2| 207 | 01852 | 383 | 37 38 37-5 
Oto 05 92 0-2392 | 22 | 23 24 23-5 
Total 968 | 959 975 | 967 


Throughout log NV has been taken as the logarithm of }(N,+N,), which 
is of course arbitrary. But we get almost equally good agreement if 


we use the expression 2(Am)? 

SAN’ 
Az being the observed number of primes in the interval AN, or the 
expression 2AN 

S,(log 


putting Av = AN /log N. 
Between 1,000,000 and 8,000,000 the agreement is almost equally 
striking as is shown in Table 1(a). Here the calculated values are 
9 
obtained from the formula Tsge 


3695.17 E 
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TABLE I (a) 
F, (eale. F, (obs.)* 
1,000,000 — 1,100,000 686-6 725 
2,000,000 — 2,100,000 625-0 644 
6,000,000 — 6,100,000 541-4 545 
7,000,000 — 7,100,000 531-5 525 
8,000,000 — 8,100,000 522-3 518 


A comparison given by Hardy and Littlewood of primes between 
100,000 and 1,000,000 is equally satisfactory. 


120 


n 
o 


2 6 0 22 26 30 38 #42 #46 S50 


AN VAY \/ 


Fic. 1 
The ordinates represent the frequencies of the various intervals 
which are plotted along the abscissa. The calculated frequencies 
are joined by straight lines; the observed frequencies are shown 
as separate dots. 


To check the formula for W, the frequencies of the intervals less 
than 100 between the primes in the region 250,000 to 255,000 are 
shown in Fig. 1 and Table 2. In actual calculations with a particular 
M it is clearly best to use the value for s,, corresponding to that 
interval rather than S,, which is the limiting value as n) > 00. The 
frequencies F,, for the various values of n are calculated by multi- 
plying the total number of intervals counted by 


p—l 
p—2 
pin, post 
n=98 
and dividing by 


n=2 p\|n,p>3 
The value of s,,, for %) = 98 is 1-475. 
* Glaisher, Messenger of Math. 8 (1879) 28-33. 


120 
; 100 
80 
40 
20 
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TABLE 2 

n F,, (calc.) F,, (obs.) | F,, (obs.)—F, (cale.) 

2 41-5 36 —5-5 6 

4 41-5 38 —3-5 6-2 

6 83 92 9 9-6 

8 41-5 41 0-5 6-4 
10 55-4 49 —6-4 7 
12 83 79 —4 8-9 
14 49-8 60 10-8 7-7 
16 41-5 37 —4:5 6-1 
18 83 82 —l 9-0 
20 55-4 51 —44 71 
22 46 46 0 6-8 
24 83 88 4 9-4 
26 46 42 —4 6-5 
28 49-9 51 1-1 71 
30 110-7 116 5:3 10-8 
32 41-5 38 —3-5 6-2 
34 44-3 51 6-7 7-1 
36 83 88 5 9-4 
38 43-9 46 2-1 6-8 
40 55-4 53 —2-4 73 
42 99-6 90 —9-6 9-5 
44 46 51 5 71 
46 43-5 48 4:5 6-9 
48 83 80 3 9-0 
50 55-4 63 7-6 8-0 
52 45-3 53 77 7:3 
54 83 80 —3 9-0 
56 49-8 47 6-8 
58 43-1 36 —T71 6 
60 110-7 108 —2-7 10-4 
62 42-9 37 —5-9 6-1 
64 41-5 39 6-2 
66 92-1 88 —41 9-4 
68 44-2 49 —5:8 7 
70 66-5 64 —2-5 8 
72 83 72 —l1l1 8-5 
74 42-7 40 —2:7 6-3 
76 43-1 41 —2-1 6-4 
78 91-8 81 —10-8 9 
80 55-4 60 4:6 7-7 
82 42-5 43 0-5 6-6 
84 99-6 97 — 2-6 9-9 
86 42-4 39 —3-4 6-2 
88 46 36 —10 6-0 
90 110-7 105 —5-7 10-2 
92 43-5 43 —0-5 6-5 
94 42-4 41 —1-4 6-4 
96 83 91 8 9-5 
98 49-8 57 7-2 7-5 


It will be observed that the mean deviation between observed 
and calculated frequencies is considerably less than VF,. 


ag 

+4 
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A better check is perhaps obtained by comparing observed and 
calculated frequencies over groups of n containing the same prime 
divisors. Table 2 (a). shows the mean values calculated and observed 
for the more numerous groups of this kind starting with the even 
values of n containing only the divisor 3, then 3 and 5, ete. 


TABLE 2 (a) 


Odd prime 


divisors of n 3 3-5 3-7 5 7 1] 0 
Number 752 329 =| «187 273 215 133 229 
F, (cale.) 83 110-7 99-6 55-4 53-7 46 41-5 
F,, (obs.) 83-5 | 1093 | 93-5 54-6 49-8 44-3 38-2 
A (in. %) 0-6 —13 | —61 —1-4 —6-5 —3-7 —8-0 
F7? (in. %) 3-6 5-5 | 7:3 6-5 6-8 8-7 6-6 


It is clear that, as the numbers go up, the percentage deviation 
goes down as it should do. The fact that the last column, which gives 
the frequency of the n’s which are powers of 2 and contain no prime 
factors, is rather lower than calculated does not seem significant. Thus 
the intervals 2 and 4 which seem abnormally low in this count between 
250,000 and 255,000, i.e. 36 and 38 instead of the calculated 41-5, are 48 
and 46 in the region 245,000 to 250,000 and average 42-1 and 42-7 re- 
spectively for intervals of 5,000 over the range from 220,000 to 255,000. 

Similar reasoning can be extended to more complicated cases. 
How for instance does the probability of P+-3p being prime compare 
with the probability of, say, P+-8 being prime if we know that both 
P and P+ 2 are primes? If P and P-+-2 are primes, obviously P+-4 
and any number of the form P+-6v+-4 must be divisible by 3. Just 
as we have hitherto omitted odd values of n, which would auto- 
matically make P+-n even, we may in this section confine ourselves 
to a sequence of values of n of the form 6v and 6v-++-2. Clearly there 
will be } as many terms in this sequence as there are whole numbers. 
But, as we need not bother about prime divisors of n which are less 
than 5, since numbers divisible by 2 and 3 have now been accounted 
for, it is clear that on the average one term in p in our sequence 6y or 
6v-+-2, must be divisible by any prime p greater than 3. Between P 
and P+3p, therefore, neither of which can be divisible by p, there 
are p—1l terms, one of which must be divisible by p. If we knew 
nothing about P+ 2, we should say that the chance of its being P+8 
was 1/(p—2). But, if we know that P-+-2 is prime, this number must 
be ruled out as a possible multiple of p, and the number divisible by p 
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must be one of the p—3 others. In other words, the chance of P+-8 
not being divisible by p will be (p—3)/(p—1) as compared to 1 in 
the case of P+-3p, which is certainly prime relative to p. Generaliz- 
ing we may therefore say that the chance of P-+-n being prime if we 
know that P and P-+-2 are primes is 

( 
times as great as the chance of P+-n being prime when m has no 
prime divisors greater than 3. 

The appearance of p—3 in the denominator in this formula of 
course indicates that, if P and P-+-2 are primes, the question whether 
a number is divisible by 3 is ipso facto decided. 

Let us count the number of prime triplets in the neighbourhood of 
N of the form P, P+2, and P-+-n, where n is any number smaller 
than m, and my is large compared with log N. According to the prime 
number theorem there will be about ,/log N primes between 0 and 
nm, and therefore n,/log N such triplets for each pair, P and P-+-2. 
But, as we have seen, the probability varies according to the number 
of different prime divisors of n. The fraction having any given value 
n will therefore be 2 


p\n, p>5 


> TT (+523) 


n=6 pln, p>5 


The denominator of course must be equal to the mean value of 


I] 


for all values of n between 4 and n, multiplied by the number of 
terms which is 4% . 
On the same lines as the above it is easy to show that 8,, the mean 


value of | | hie. +5"5) converges to 
1 


1-38 


6 ... for large values of no. 


= 


which is abou 


= 
‘ 
| 
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As has been shown, the number of pairs in an interval AN is 
2 AN 
S, log NV log N° 

Each forms about ”,/log N prime triplets in which n in the term P-+-n 
is less than m). The probability of finding one of them having the form 
P, P+2, P--n is, as we have seen, 

S30 p—3 
Hence, for large values of m), the number of triplets of the form 
P, P+2, P-+-n in the interval AN should be 


6 2 
14 


p>5 


For n = 6 we find simply 
6AN 2-86...AN 


Since all our arguments apply equally to positive and negative 
values of n, it is clear that the frequency of the triplets P, P+-2, P+6 
must be the same as that of the triplets P, P+-4, P+6. 
In Tables 3 and 3(a), in which the sum of the frequencies of the 


TABLE 3 

F, (cale.) F,, (obs.) 
22-9 20 
10-0 7 
14:7 19 
31-7 31 
40-6 29 
119-3 113 
90-110 75:0 83 
220-255 105-5 109 
419-7 411 

TABLE 3 (a) 

AN.10-5 F, (calc.) F (obs.) 
339-8 333 
298-4 312 
275:1 312 
259-3 248 
247-6 244 
238-4 248 
230-9 220 
224-6 200 
219-2 215 


2333-3 2332 


| 
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two types of triplets is used, observed and calculated values are 
compared. Our constant 
p*(p—2)(p—3) 3p?—1 
I] (p—1)\(p—2) I] 

is the same as that found by Hardy and Littlewood from whose 
paper (loc. cit.) the number of triplets shown in Table 3 (a) is taken. 

The numbers in the two groups are reasonably nearly equal. In 
Table 3 there are 202 of the type P, P+2, P-+-6 as against 206 of 
the type P, P+4, P+6. In Table 3(a) the corresponding numbers 
are 1,138 and 1,194. 

A similar though more complicated argument would give the 
formula for the frequency of quadruplets.* 


* Though not strictly relevant it is perhaps interesting to point out that 
the number of primes of the form 4m?-+-1 can be derived by similar methods. 
The number of integers of this form in an interval AN is }AN/VN; the proba- 
bility of an odd integer being prime is 2/log N. Hence the number of primes 
of the form 4m?-+1 less than n should be proportional to 4Livn. To take 
account of the fact that 2 out of every p terms in the series 4m*+ 1 are divisible 
by p we must multiply this expression by 


p, being primes of the form 4v+1. To allow for the fact that primes of the 
form p, = 4v+3 cannot divide 4m?-+ 1 we must multiply by 


If we do this, we arrive at the same formula as Hardy and Littlewood, namely 


IT] 


Their formula for the frequency of prime pairs of the form 4m?+1 and 
4m*+3 may be found on similar lines. 

The number of representations of any even number as the sum of two 
primes (Goldbach’s Theorem) may also be derived by probability arguments. 
If the number of different prime divisors of n is small compared with the 
number of primes less than n, the probability that n—p, will be a prime, n 
being a number re N and N+ AN, will tend as above to 


TES | 
p—2 
p\n 


if we average for all values of N less than n. The number of primes P, less than 
n is Lin. Hence the probable number of occasions on which P,+P, = n 
p—2 +5 (log n)? 


as found by Hardy and Littlewood and by as 


=; 
Pi~ 
D323 | 
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The results derived from the assumption that the distribution of 
the primes may be treated as random are so well borne out that it 
seems worth while to carry our considerations a stage further and to 
see whether it is possible, albeit on very rough and ready lines, to 
calculate the frequency with which intervals between neighbouring 
primes occur. 

Between P and P-+-n there are }(n—2) odd numbers. Normally, 
therefore, if we call the chance of P+-2v not being a prime ay,, v 
being any whole number, the chance of no prime intervening between 
P and P-+n, should be 


The probability w,,, therefore, that we shall find two neighbouring 
primes separated by an interval n should be equal to W,,, the inherent 


probability of this interval being found between two primes at all, 


n? 


multiplied by 


v=1 


which is the probability that all of the intervening numbers are 
composite. 

It remains to find an expression for a,,. The first approximation is 
clearly given us by the prime-number theorem. If the chance of any 
number being a prime is 1/log NV, the chance of an odd number being 
a prime is 2/log NV, and the chance of an odd number not being a 
prime consequently 1—2/logN. Therefore we may expect to get 
something like the right result by putting 

2 
A, = — 
log A 

* Tt should be pointed out that this rough and ready approach is not 
strictly consistent with the formula derived above for the frequency of 
triplets of primes of the form P, P+2, P+6 or P, P+4, P+6. If we take 
no account of intruding primes, there should be twice as many pairs of the 
form P-+6 as of the form P+2 since W, = 2W,. Hence, if the probability of 


/ 2 
P +6 being a neighbour of P is only 217,( 1— icy) it would follow that the 
g 


frequency of triplets should be 
logN “logN 8S, (log 
instead of 12 x 
(log N)* 


as derived more rigorously above,—a difference of about 8 per cent. 


Vv PALL < 
| | 
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This assumption gives us the same value for the frequency of prime 


2AN 


pairs as we derived above, namely Slog.) If we designate the 


probability of intervals between neighbouring primes by small letters, 
clearly W, = w,. Further, if «,, = 1—2/log N, 


\Hn-2) 
wv, = ( + 
log V p—2 
Since P must be followed eventually by some prime, s w, = 1. 
2 


1 
value S,, as we may by a well-known theorem,* we find 


2 2 \3 S,w 


which gives us for the probability of a prime P being followed by 
another P+ 2, 


Substituting for the fluctuating expression I] (1 + 5) its mean 


9 


= Blog N’ 


as above. Obviously the crude expression 


2 4(n—2) 1 
p\n, prs 


gives us also the right value for the mean value % of n. For, if we 
again take the average value of 


it is clear that 2 nn = logN. 
> w, 

In Table 1, F, the frequency of the interval 4 was given as well as 
F, and indeed was averaged with it in the last column. The reason 
is, of course, as we saw when considering prime triplets, that w, must 
be equal to w,. For it is just as impossible that a prime should inter- 
vene between P and P-+-4 as that a prime should intervene between 
P and P+2. No prime can intervene between the two odd primes 
P and P+2 because the only intervening number is even—i.e. 


* Hardy and Littlewood, Proc. London Math. Soc. 13 (1913), 174. 


eee 

as 
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divisible by 2. No prime can intervene between P and P-+-4 because 
the only odd number between P and P+-4 is P+2, and P+2 must 


99 


22-26 
Fia. 2 (220,000 to 255,000) 


In these diagrams (Figs. 2-4) the ordinates represent the frequencies of the 
various intervals which are plotted along the abscissae. The calculated fre- 
quencies are joined by straight lines; the observed frequencies are shown as 
separate dots. 


500 


Fic. 3 (10,000 to 30,000) 


be divisible by 3 if P and P-+4 are primes. In Table 4 and Figs. 2, 
3, and 4 observed and calculated frequencies of the intervals between 
neighbouring primes are compared. f, (cale.) and f, (calc.) are put 
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SOF 
30 
20+ 
10+ 
Fic. 4 (500 to 2,000) 
2A 
equal to TineW¥ as above; all the other terms are calculated accord- 
2 108 + 


ing to the formula 

i+) 

without any arbitrary constant. They show surprisingly good agree- 
ment between counts and calculation in the intervals 220,000 to 
255,000, 90,000 to 110,000, 10,000 to 30,000, 2,000 to 10,000, and 500 
to 2,000. This is the more remarkable since our formula for W, 
ignored all second-order terms. In w,, we have introduced a term in 
1/(log NV)? although we had neglected the possibility of such a term 
in W,, which would of course appear in w,,. That the discrepancies are 
not larger seems to indicate that the second-order terms in W, are 
relatively unimportant. 

Although the observed and calculated values in these tables seldom 
differ by more than the expected statistical deviation, careful scrutiny 
does reveal a tendency to certain systematic discrepancies. In 
general the observed values seem to fall off rather more rapidly 
towards the end than the calculated. But there also seems to be 
a more specialized trend. Had we not realized that w, must be equal 
to w, and left the crude value . 


we should have found observation systematically exceeding prediction 


3 
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TABLE 4 
220,000 90,000 10,000 2,000 500 
to to to to to 
255,000 110,000 30,000 10,000 2,000 
2 299-5 | 295 | 199 198 267 262 140°5 | 144 38-5 37 
4 299-5 299 199 | 206 267 267 140-5 141 38-3 38 
6 502 485 329 | 312 | 426 426 218 219 55-2 55 
8 210 197 135-8 | 135 170 168 84-5 74 19-8 19 
10 235 271 149-6 | 174 181-3 | 167 76-2 91 19-2 22 
12 296 292 | 185-5 | 157 | 216-8 | 220] 101-3 84] 20-7 16 
14 148-5 153 91-8 97 94-1 103 47-2 43 8-9 7 
16 103-8 | 103 63-4 64 69 75 30-5 32 5:3 2 
18 174-2 159 104-7 103 110 110 47-3 34 7-7 6 
20 97-2 87 57-7 48 | 58:8 48 24-5 13 3:7 2 
22 68 79 39-7 43 39 42 15-8 14 2-2 2 
24 102-5 101 59-1 7 6-1 30 22 14 3 1 
26 47 50 26-6 24 24-5 12 9-3 3 1 
28 43°3 43 24-2 20 21-5 22 8-0 5 
30 80-7 57 44-3 32 38 27 13-6 11 
32 25-4 17 13-8 7 11-4 3 4-0 1 
34 22-7 14 12-1 13 9-7 9 3°3 2 
36 35°7 25 18-8 15 14:5 ll 5-2 1 
38 15:8 14 8-2 5 6:1 0 2-0 
40 16-7 | 14 8-6 3 6-6 2 
42 25-1 | 19] 12-7 5 
44 98 | 7 49 3 
46 77) 3 4-0 5 
48 9 6-8 2 
50 9-3 3 ] 
52 4:7 0 2-2 | ] 
54 73 | 3 3-4 3 
56 3°5 | 1 1-7 
58 2-6 | 2 | 
60 5:7 5 1 
62 1-9 0 


at n = 4. The other main systematic deviations, namely a tendency 
to an excessive number of intervals at n = 10, n = 16, and n = 22 
are due to the same cause operating less obviously. 

If we take divisibility by 3 into account, the number of places 
where primes can intervene is clearly substantially reduced. In the 
following table the second line shows /,, the number of places in 
which primes might intrude having regard only to divisibility by 2, 
the third line the number if both 2 and 3 are considered. 

n= 24 6 8 10 12 14 16 18 20 22 24 


L= 0 1 


l= 0 


23 4 
12 2 


5 6 7 
3 4 4 
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9 #10 


5 6 6 


11 
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Clearly the number of possible intruders instead of rising steadily as 
in the second line halts always when n = 6v-+-4, v being any whole 
number. It is easy to take this into account when /, = 0. The 
probability of P+-2 being composite is 1 if P and P+4 are primes. 
But the answer is not so simple when /, halts at 2 instead of proceed- 
ing to 3. For the term 1—2/log N, which gives the probability that an 
odd number will be composite, must be modified. We cannot just 
substitute 3/log N for 2/log N. True, this represents the probability 
of an odd number not divisible by 3 being prime. But this is not 
what we have got. We are dealing with the probability of a number 
between P and P-+-n being prime, knowing P and P-+-n to be prime. 
And to try to take this into account quantitatively is very com- 
plicated, as we have seen even in the comparatively straightforward 
case of prime triplets of the form P, P+2, P+6. Qualitatively, of 
course, it is easy to see what the result will be. Instead of the 
uniform exponential decrease there will be at these values of n a 
slowing up, leading to a slight excess of intervals between primes 
when n = 6v-+-4 over the number calculated by the simple expres- 
sion. And this is what we observe. 

This note does not claim to be more than a first step towards 
calculating the distribution of the intervals between primes by 
statistical methods. Clearly various other more intricate formulae 
could be derived by somewhat more elaborate processes and very 
likely better agreement could be achieved between calculation and 
observation. But even this crude approach gives results so closely 
in accord with observation that it would be difficult without an 
immense amount of statistical material to establish whether a 
genuine improvement had been made. As always in considerations 
of this kind it is questionable where the line should be drawn, i.e. 
where numbers are large enough to justify simple statistics. 

It is, of course, possible that we have been fortunate in the 
selection of the relatively small sample it has been feasible to 
examine and that closer and more detailed study will reveal dis- 
crepancies not obvious in this preliminary investigation. But that 
these methods are capable of yielding interesting results seems to be 


proved. 

Above all it seems remarkable that the formulae such as those for 
the frequency of pairs and triplets derived by much ‘deeper’ methods 
by Hardy and Littlewood based on an extension of Riemann’s 
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hypothesis, should arise quite naturally in such a very simple way. It 
would appear interesting to examine whether and in what manner the 
Riemann hypothesis and Hardy and Littlewood’s extension thereof 
are linked with the hypothesis that the distribution of the primes 
may be treated as random. 

In conclusion I must express my thanks to Mr. J. Harvey for the 
care and skill with which he has carried out the counts which have 
enabled the comparisons to be made. 


Nore ON THE PREvIoUS WorRK 


Hardy and Littlewood use their celebrated ‘Partitio Numerorum’ 
method to find likely asymptotic formulae for the number of pairs of 
primes (P, P+2n) less than a large nm and similar formulae for triplets, 
quadruplets, and soon. Their method is entirely distinct from that used 
in this paper and ‘deeper’; the assumptions they make are very different 
in form, but fundamentally not dissimilar. Staeckel, on the other hand, 
uses probability arguments about as ‘elementary’ as in this note. His 
papers are long and, like Hardy and Littlewood’s, deal also with Goldbach’s 
Theorem ; it is difficult to pick out exactly those parts which correspond 
to the work set out above, but the methods here used are much shorter 
than his. 

Results equivalent to those found above for W,, were obtained by Hardy 
and Littlewood. They do not deal with the problem of w,,, though Staeckel 
gives certain preliminary considerations. In the triplets problem Hardy 
and Littlewood find a formula identical with that found above. Staeckel 
seems to have a result similar to theirs. 

Hardy and Littlewood remark that a better comparison with numerical 


n 


5 in place of 


data is always given by using | lees} il - 
og x)* ogn 


; in the formula 


for the numb.: of prime-pairs less than n, and similarly in the other 
problems. The results derived by the methods used in this note appear 
naturally in terms of ‘frequency’, but they correspond, of course, to the 
integral form. It is presumably for this reason that the neglected second- 
order term in W, turns out to be so small. A further advantage of the 
frequency form is that it enables numerical comparisons to be made over 
any interval; we do not have to start from 0 with our count. 

So far no impression has been made on these problems by rigorous 
analysis. One of the best results so far proved in this direction appears to 
be that due to Buchstab, Rec. Math. (Mat. Sbornik), N.S. 10 (52), (1942), 
87-91, who proves that, for any A, there exist infinitely many primes p such 
that each prime factor of P+2 exceeds (log p)*. 


ELECTRICAL NOTES 
By F. B. PIDDUCK (Ozford) 
XII. ALTERNATING CURRENTS IN NETWORKS} 
[Received 27 October 1944] 


THE mutual impedance of two open circuits at any frequency has 
been found by Murray.{ Let P be a point on a circuit distinguished 
by the upper affix uz at a distance p from a fixed point on it, measured 
along the arc, and let the ends correspond to p = p, and p = py. 
Let Q similarly be a point on a circuit v distant g from a fixed point 
on it, and let PQ = r. Let currents of frequency w/2z flow in the 
circuits, and let the current at P be I#(p)expiwt, where J#(p) is in 
general complex and the real part of the product is taken for a 
physical interpretation. Let J(p) when split up into harmonics be 


of the form - 
>, 


Murray found that the mth and nth harmonic currents in the circuits 
p and v have a mutual impedance given by 
Ds 

where we have used rationalized metre-kilogram-second units and 
written k for 27/A and y for the angle between the tangents at Pand Q. 
To include both open and closed circuits, we have used the conjugate 
complex J#*(p) instead of J#(p) itself, and in this notation the ten- 
sion applied to the circuit pu is 


Py 
Uy, = | dp, 
Py 


where E#(p)exp iwt is the component of the electric intensity at P 
parallel to the tangent, in so far as it is applied from outside. 

By considering Z##, we find the self-impedance of a single circuit 
of thin wire, considered by Murray, and also extend the formula to 
networks. The double integral is not strictly a double line integral, 

+ This note is a continuation of Note IV, Quart. J. of Math. (Oxford), 


2 (1931), 174. 
t F. H. Murray, American J. of Math. 53 (1931), 873. 
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m 


distributed over the surfaces of the circuits. In calculating Z4# we 


mn 


but a double surface integral, since the currents J/#(p), I},(q) are 


have to take the surface integrals twice over the same surface. Con- 
sider the contribution of the parts of the integral where r is small. 
It will not affect the argument much to take the current to be dis- 
tributed uniformly over the cross-section, supposed circular. When 
r is small, exp(—ikr)/ikr is nearly equal to 1/ikr, and the contribution 
is [ | dSdS' /ikr, where a constant has been dropped and dS, dS’ are 
elements of surface. This is proportional to the potential energy of 
a uniform charge on a cylinder, which is equal to the mutual potential 
energy of two line charges, one in the line of centres and the other 
parallel to it in the surface of the wire. Thus Z##, for a single wire 
is Z#”, for two such parallel geometrical lines, or filaments. Similar 
considerations apply to a network. Let the wires » and v have in 
common a piece of wire AB. Then in calculating Z4”, we must take 
the arc of one circuit there to be the line of centres and the arc of 
the other a parallel filament in the surface of the wire. 

This theorem is more general than my former theorem in that the 
current need not be the same at all points, but less general in being 
restricted to currents lying wholly in the surface. It applies equally 
to closed or open circuits, except that J#(p) for a closed circuit need 
not vanish at any particular point of the wire. It is common to 
excite an antenna by attaching a feeder to two points on either side 
of the centre, and I have shown by these formulae that the best 


position of the points can be found approximately by theory. 
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